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[1] Examenes con DERIVE , Jose Luis Llorens Fuster 

DERISOFT (FAX: 907 330 371), 214 pages, ISBN 84-921694-1-9 

|2] DERIVE: Aplicaciones matematicas para PC, Carrillo & Llamas 
RA-MA, Madrid, 1994 

|3] Como hacer Matematicas con DERIVE, Castro Chadid 
Reverte Colombiana, 1993 

Aplicaciones de DERIVE: Analisis Matematico-I (Calculo), J.L. Llorens Fuster 
Servicio de Publicaciones de la Universidad Politecnica de Valencia, ref. 93735. Valencia, 1993 

Aplicaciones de DERIVE: Algebra Lineal (Fundamentos), J.L. Llorens Fuster 
Public, de la Universidad Politecnica de Valencia, ref. 94753. Valencia, 1994 

Introduccion al uso de DERIVE: Aplicaciones al Algebra Lineal y Calculo Infinitesimal, 

J.L. Llorens Fuster, Public, de la Universidad Politecnica de Valencia, ref. 95638. Valencia, 1995 

[7] Calculo Matematico con DERIVE para PC, Paulorrogan & Perez, 

RA-MA, Madrid, 1994 

[9] Analysis mit DERIVE, Lehr- und Arbeitsbuch fur die Jahrgangsstufe 11/1, Hans Jurgen Kayser 
Diimmler Verlag, 1996, 157 Seiten, ISBN 3-427-4523 1-X; Diimmlerbuch 4523 

Ich erhielt mehrere Anfragen, da das Buch von Larry Gilligan "Mastering the TI-92" fur den deutschen 
Sprachraum offensichtlich schwer erhaltlich sein soil. Die ISBN ist 0-9626661-9-X, erschienen bei GILMAR 
Publishing, P.O.Box 6376, Cincinnati, OH 45206, USA. Sie konnen das Buch aber noch bequemer erwerben bei 

bk-teachware, FAX ++43 (0)7236 6065 80. 
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Exchange for DERIVE - Teaching 
materials in the DNL. 

The wHeeC Has not to 6e invented twice. 


Borse fiir DERIVE Unterrichts- 
materialien im DNL. 

(Das (Rad muss nicfit zweimal erf unden werden. 


I can offer: 

Binomial Theorem, GCD & LCM, System of coordinates, Modelling Word problems with 
DERIVE (all in English and German/, SET. EXE, MENGE.EXE (a program to teach and to 
practise operations with sets including VENN diagrams, in English and Gennan). 

I have produced a paper „Einfuhrung in die Matrizenrechnung mit DERIVE “ - in German. If 
you are interested I would send you the paper on a diskette in MSWord6-format. 



This is DfW 

DERIVE for WINDOWS. 


The first screen shot using DfW, de- 
monstrating an ALGEBRA- and a 
PLOT window. 


Cm;;;;: -0.9189?, -1.3S?9 


Center: 0. 


Scale: 1:1 


One of the new features: You are 
able to add annotations to your plots. 













Liebe DUG Mitglieder, 

Es war wirklich schdn, so viele von Ihnen in Bonn zu tref- 
fen oder personlich kennen zu lernen. Fast alle, die zum 
ersten Mai erne DERIVE Konferenz besuchten, waren 
einmal mehr begeistert von der familiaren Atmosphare, 
vom hohen Niveau und vom ungebrochenen Enthusiasmus 
der Teilnehmer. Sie haben es genossen, mit David Stoute- 
myer, mit Bert Waits, um nur zwei zu nennen, und vielen 
anderen, die sie aus dent DNL kennen, zu sprechen und zu 
diskutieren. Ich wurde einige Male angesprochen: "Josef, 
zeige mir bitte den ..., ich mochte doch wirklich wissen, wer 
das ist!" 

Nicht nur J.Wiesenbauer nahm einige 
Ideen von Bonn nach Hause. Auch Josef 
Lechner und Eugenio Roanes haben ihre 
Kontakte vertieft und sind nun dabei, 
gemeinsam, LOGO in DERIVE zu imple- 
mentieren. Leider konnte ich bei so vielen 
Vortragen nicht, die mich interessiert 
hdtten nicht dabei sein, aber dafiir freue 
ich mich auf den Tagungsband. 

Abschliefiend mochte ich namens der DUG 
den beiden Organisatoren, Barbel Barzel 
und Leo Klingen fur Ihre grofie Miihe 
herzlich danken. 

Vor Bonn hatte ich die Gelegenheit, an einer T 3 -Ausbildung 
- Teachers Training for Technology - in Columbus, Ohio, 
teilzunehmen. Diese, von B. Waits und F.Demana gegriinde- 
te Initiative soil nun unter anderem auch auf Europa aus- 
gedehnt werden. Bei dieser Gelegenheit darf ich Sie ersu- 
chen, TI-92 Beitrage an Bert Waits oder an mich zu schi- 
cken. (Vergessen Sie aber bitte nicht das klassische 
DERIVE!) 

Thomas Weth kommt mit seinem Kurvenlexikon wieder im 
nachsten DNL zu Wort. Gerade bei seinen Beitragen lasst 
sich eine sinnvolle Verbindung zwischen einem CAS und 
der interaktiven Geometrie sowohl auf dent PC als auch 
auf dem TI-92 herstellen. 

Ich mochte Sie nochmals darauf hinweisen, dass Sie die 
.MTH-files zu diesem DNL von SWHH's web page herunter 
laden konnen. Danke, Al Rich. 

Vom renommierten Verlag BIRKHAUSER AG wurde ich 
gebeten, fur unsere deutschsprachigen Mitglieder eine 
Information iiber die Fachzeitschrift "Elemente der Ma- 
thematik" beizulegen. Diesem Wunsch komme ich gerne 
nach. 

Bis zum nachsten DNL - der „ Weihnachtsnummer" mit 
einigen Zugaben - griifit Sie herzlichst 

Josef Bohm 
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Dear DUG Members, 

It was really nice to meet so many of you in Bonn and 
to make your acquaintance. Most of those who at- 
tended a DERIVE Conference the first time were 
fascinated from the familiar atmosphere, the high level 
and the unbroken enthusiasm of the participants. They 
enjoyed talking, discussing and joking with David 
Stoutemyer and Bert Waits, to mention two people, 
and with many others who they were knowing from 
the DNL. It happened several times that I was asked 

"Oh, Josef, please show me , / really would like to 

know who he (she) is!" 


Not only Johann Wiesenbauer took 
some ideas back home. Josef Lechner 
and Eugenio Roanes for example, 
have renewed their contacts and are 
now busy to implement LOGO in 
DERIVE. Unfortunately I was not able 
to attend all the lectures I was inter- 
ested in. But I am looking forward to 
studying the Proceedings. 

Finally I want to express many thanks 
on behalf of the DUG to the organizers, 
Barbel Barzel and Leo Klingen for their 
tremendous efforts to make our stay in 
Bonn as pleasant as possible. 


The week before Bonn I had the possibility to attend a 
T 3 - Conference - Teachers Training for Technology - 
in Columbus, Ohio. T 3 , founded by B. Waits and 
F.Demana will be established now in Europe, too. On 
this occasion I’d like to invite you to submit TI-92 pa- 
pers to B. Waits, B.Kutzler or to me. (But don’t forget 
the classic DERIVE!) 

Thomas Weth will continue his "Lexicon of Curves" in 
the next issue. Especially his articles are excellent 
examples how to link a CAS with interactive geometry 
which can be realized now on a PC as well as on the 
TI-92. 

I want to point out once more that you can download 
the .MTH-files belonging to this issue from SWHH's 
web page. Many thanks, Al Rich. 

Until the next DNL - the "Christmas Issue" with some 
extras - 


Sincerely yours 


Josef Bohm 
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T SI II f| I* SI III K 
with nr . hive 

(in one of the next ACDCs) 


Nur fiir Eingeweihte: Habt Ihr das Mdppchen gefunden? 
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The DER1VE-NEWSLETTER is the Bulle- 
tin of the DERIVE User Group. It is pub- 
lished at least four times a year with con- 
tents of 40 pages minimum. The goals of 
the DNL are to enable the exchange of ex- 
periences made with DERIVE as well as to 
create a group to discuss the possibilities 
of new methodical and didactical manners 
in teaching mathematics. 

We include now a section dealing with the 
use of the TI-92. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/FAX: 43-(0)660 3136365 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE Newsletter will be. 


Preview: Contributions for the next issues 

Graphic Integration, Linear Programming, Bohm, AUT 

LOGO in DERIVE, Lechner & Roanes Lozano, AUT & ESP 

3D-Geometry, Reichel, AUT 

Parallel- and Central Projection, Bohm, AUT 

Algebra at A-Level, Goldstein, UK 

Tilgung fremderregter Schwingungen, Klingen, GER 

A Utility file for complex dynamic systems, Lechner, AUT 

Examples for Statistics, Roeloffs, NL 

Linear Mappings and Computer Graphics, Kummel, GER 

Solving Word problems (Textaufgaben) with DERIVE, Bohm, AUT 

DERIVE and ACROSPIN, Schorn & Bohm, AUT/GER 

Line Searching with DERIVE, Collie, UK 

About the "Cesaro Glove-Osculant", Halprin, AUS 

3rd-order Differential Equations, Guitierez, ESP 

Tangrams with DERIVE, Poblacion, ESP 

The TI-92 Section, Waits a.o. 
and 

Setif, FRA; Vermeylen, Belgium; Leinbach, USA; Halprin, AUS 
Weth, GER; Wiesenbauer, AUT; Keunecke, GER; Aue, GER; 
Weller; GER, Stahl, USA; Mitic, UK; Koth, AUT; and 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 









We have recently met at the DERIVE Conference in Bonn. 1 liked the atmosphere very much and I have found a 
few colleagues teaching the same subjects as I do. My favourite interest is financial math. 1 have recently 
visited the ICME at Sevilla and last year I was at the CIAEM Conference at Berlin and there was no lecture on 
financial math. The subject seems too practical for mathematicians and too difficult for economists. So my ques- 
tion is: could we start a user group of financial math teachers, whatever the platform? Personally, 1 have tried to 
do this in Holland - with Hewlett-Packard financial calculators - with no success. People are just not interested 
in any machine in the classroom. By the way; there are no graphic calculators and there is no DERIVE in Dutch 
secondary education, and there will not be anything like that on the central examinations for the next four years. 


DNL: . I think we could establish a "Financial Math Working Group" within the DUG if there are other 
interested members, too. Or we start the Group as a pair of enthusiasts. I would like to ask you for a 
"pilot contribution" to make clear some topics for this group. As I am a teacher on a "Handelsa- 
kademie" (- "College for Business Administration") I am very interested in your idea and I am looking 
forward to hearing from you in the near future. 

Jan Vermeylen, Kapellen, Belgium 

Have you got any experience with the link you can make between DERIVE and the programming language 
C++? Where could 1 find more information or examples? 

DNL: Are there any members with experience in combining DERIVE results and C++-applications. I 
used very successfully the linkage to BASIC programming ACD.EXE to convert DERIVE results into 
ACD-files for ACROSPIN. I’ll present this program in the next DNL. Any other applications would be 
appreciated. 

George Feeman, A1 Ain, United Arab Emirates 

We have come across this bug in Ver.3.00 and 3.10: 

2 2 6 


3 + x 

X 

Simplifies to 

x(x + 3) 

2 

3 

Simplifies to 

x + 7 

and so on, but 

(x + l)(x + 3) 

3 + x 

x + 1 

2 

3 + x 

3 

X 

DOES NOT Simplify - the output is the same expression. 


Can you tell me why? Why this exception to this simple algorithm? 

DNL: At the occasion of the Bonn DERIVE Conference I discussed this "bug" simplifying the sum of 
two simple rational expressions with David Stoutemyer. By the way: that’s not a "bug" of only DERIVE 
3.x versions. DERIVE 2.58 returns the same results. That’s a summary of David’s answer: 

Computer Algebra is among others the art of recognizing and using patterns. DERIVE always looks 
for patterns to apply any simplification. (Josef’s guestion: "What is the definition of simple???") 

2 2 

In DERIVE recognizes the same numerators and then tries to apply a rule. 

3 + x x 


3 + x x + 1 


- DERIVE discovers two denominators of the form (x + constant) and then ap- 
plies the rule to add two fractions using common denominators, but 


2 3 3 „ _] 3 

DERIVE cannot find any pattern: — is treated rather as 3 x than as 

3+xx x x+c 


Maybe you will not find this explanation very satisfactory. Nevertheless it is always interesting to look 
behind the curtain of a CAS. 








Helmut Wunderling, Berlin, Germany 

Anbei Bildschirmausdmcke von DERIVE. Die obe- 
ren Bilder zeigen, wie unter SoftWindows auf einem 
Powermac die Aufbauroutinen des Bildschirms nicht 
vollstandig funktionieren. Das hat zur Folge, dass bei 
Ausdriicken, die wechselnd ein- und mehrzeilig sind, 
die Lesbarkeit dahin ist (vgl. Ausschnitt). 

Der untere Teil zeigt, dass Version 3.12 die Sonder- 
zeichen fur Integrate, Wurzeln, Summen usw. nicht 
bereitstellt und zwar auch auf PC (z.B. Hewlett Pa- 
ckard) und das, obwohl die Version 2.59 z.B. das nun 
wieder anstandslos tut. 

Ich brauchte aber fur Veroffentlichungen gute Bilder. 
Was kann ich tun?? 


Attached you can find screen dumps of DERIVE. 
The upper pictures (full screen and part of the 
screen) show that using SoftWindows on a 
Powermac the building up routines of the screen 
will not work properly. From this follows that 
expressions with alternately one and more lines 
are not readable. The dump on the bottom 
shows that version 3.12 does not provide the 
special characters for sum, integral, roots etc, 
and that even on a PC (Hewlett Packard) al- 
though I don’t have any problems using 2.59. 
For a publication I need nice screen dumps. 
What shall I do??? 


n nathenatLea] flas-fesiant 


Vers Lon 3,12: 


Copyright CC1 1980 through 1096 Silt 
rehouse, Itk 

3h60 KiJfllpe fluenue, Suite 301 
hortolu In h llaua i j , 96S16-323fc, J Sfl 

lease da rmt. nakc illegal cop I as i>F ttEHlUEl Tli La software is not sharrajnre or 
freeware. U U TOt to be put I Lstiei on bulletin hoards at distributed bij nnu 
■ther r>a aTis ul th™,it or It ten pern i as inn Frnm So Ft Uarelwuse. Coe. 

or technical support. or If you kimu Of ary. perscni or caPOp<ain|L^ Is-Etr JlmtLng 
£RJir£ 4 S shareware nr freeware, 0 lease write -us at the eboUC Or send a 

{■?* to C0O81 T3S-11B5 . 

Tress W for help [jj 



Ci? ? f-H- tur 3^:^ srjZ i\: * 

■ tfr 2 : ■ 


^ r 3 2 1 

- 5 * +5 

HI SQST — A* 

S *-3 

i 

BZ: BTif (4*1 Aa> Alt) 



I h^rpi:; • iZL 
unfit tc t i fid - 
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In the following you can see how to transfer this utility to TI-NspireCAS. The rank function is not 
implemented in the base system but it is available via the linalgcas- library. The screenshots below 
should be informative enough that you can reproduce the procedure on your own. 


Catalog 

Jl Math Operators 
Unit Conversions 
[Qj Libraries 

Except for current problem, items listed below 
are in the MyLib folder within C;\Dokumente und 
Einstellungen\User\Eigene Dateien\TI-Nspire\ 

§§ Refresh Libraries 

UC4J41CIII 

diagonalization 

dn 

eigenvals 

eigenvects 

expmat 

gausstep 

help 

inversestep 
kemelbasis 
kemelvectors 
pwrm at 
rank 

simultstep 


* 


□ ^ Wizards On 


rank(Mat): rank of matrix 


c/7/erewc;t7s(v):=list^mat(seq(|v[l,^:+l]-v[l ) A:]| ( A;l ) dim(v)[2]-l)) 

Done 

1 

q larga r( v) : =au gment(au gment( [o],v),[v[l,d im( v) [ 2 ] ] ] ) 

Done 1 

dim _jordan(v)=dferencias[diferencias[aIargar[v))) 

Done 1 

1 dimnu1os[a,z]\= list>mat 

fseq(dim(<7)[ ll-fl^^SWS^^I(g-Z'identitv(dim(fl)l'll)y ? Ip.l.dimMfll)) 







Done 1 

1 bloq_jordan[a,z ): 

:=dim _Jordan[dimnulos(a,z)) 

Done 

| © Examples: 


2 0 0 0 


2 0 0 0 


2 0 0 0 


mJ\= 

12 0 0 

\m2\= 

12 0 0 


12 0 0 



0 12 0 


0 10 0 


0 10 0 



0 0 0 2 


.0 0 0 5 


0 0 0 5 


1 bloq Jordan{m 1,2 ) 

'lOlo' 


1 1 eigVl(/w2) 


{ 0 . 00000000000 , 2 . 00000000000 , 2 . 00000000000 , 5 . 00000000000 } 


1 1 Iinalgcas\eigenvals(m2) 


{ 0 , 2 , 5 } 


1 bloq_jordan[m2,'. 

2) 


[ 0100 ] 

- V 

22125 | 


Document 1 x Jordan x f 


rank 1 0/1 


Define LibPub rank(m)= 

Func 

©rank(Mat): rank of matrix 
Local n,s 

s\=check_type_arg{m,"'MA1" ,o) 
If s*"ok" Then 
s 

Elself is_cp_matrix[m) Then 
Ub_msg { " too_cplx ”, false) 

Else 

n:=rref{m) 

rowDim(ff) 

(when(norm(n[i])=0,0,l)) 

! = 1 
Endlf 
EndFunc 




bloqjordan[m2,5 ) 


m3 = 


17 1 1 _ "8 

— — 0 0 — 0 — 

3 3 3 3 

4 11 „ „ -1 „ -1 

3 3 3 3 

25_ 1_ 1 ^ ^2 x 

6 3 2 3 3 

-1 1 0 4-102 

5 1 „ „ 13 -8 

3 3 3 3 

- 1 -4 — -2 5 -1 

2 2 

5 1 „ „ 1.4 

— — 00—0 — 

3 3 3 3 


[10 0 0 

17 1 1.-8 

— — 0 0 — 0 — 

3 3 3 3 

4 11 „ „ -1 -1 

3 3 3 3 

25 _ 1 _ 4 ^ -8 

6 3 2 3 3 

-1 1 0 4-102 

5 1 „ „ 13 -8 

3 3 3 3 

- 1 -4 — -2 5 -1 

2 2 

5 1 „ „ 1 „ 4 

3 3 3 3 


!inalgcas\eigenvals(m3) 

{ 3 , 4 ) 


bloq Jordan[m3, 3 ) 

[1100000 


bloq Jordan[m3, 4 ) 

[1010000 


[ 3 2 0 0 ol [ 1 3 2 1 4l 

[l 3214 

v 

12/25 
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linalgcas\eigenvals[m5 ■ m4‘ m5 1 ) 

{5} 

Iinalgcas\ceigenvals{m5 ■ m4 ■ m5 ^ 

■) < 

3+2 i, 3-2 1,5} 

bloqJordan{m5 ■ m4 ■ m5 ^,3-2' i) 


[OIOOO] 



3 

2 

0 

-1 


1 

1 

-1 

0 


1 

1 

-1 

0 

m6:= 

-1 

1 

1 

1 

\u6\= 

-1 

0 

1 

0 


-1 

0 

1 

0 


1 

1 

2 

0 


1 

0 

0 

0 


1 

0 

0 

0 


-1 

-1 

0 

2 


-1 

0 

0 

1 


-1 

0 

0 

1 


eigVl(w<5 

{ 2. 00000003722, 1 . 99999996278,2. 00000000000,2. 00000000000 } 

linalgcas\eigenvals[m6] 

{ 2 } 

bIoqJordan[m6,2 ) 

OOOl] 

u6 1 m6u6 

2 10 0l 
0 2 10 
0 0 2 1 
0 0 0 2J 


For accomplishing this contribution see the same procedures performed with the TI-92 / Voyage 200. 


(F 17 *H)T F 2 T Y F 3 t Y FH Y FE Y 1 

1 t f— Command |U i eu lExecute IF i nd... 1 1 


fF 17 ^H)Y F£t Y F 3 t Y FI Y FE Y i 

1 t f— ICommand |U i eu lExecute IF i nd... I 

:0 aiferenciae 



C: 1 i£tFmat(seq(abs(u [ 1 , k +1 ]-u [ 1 , k 1 5 ,fc, 1 , 
dimCul [ 2 ]-lU) 4 diff<vJ 


:0 dimnulos 

C:: 1 istFmat(seq(dim(al [ 1 ]-rank( (a-z*ident 

: 0 a 1 ar gar 

C : augment (augment ( [ [ 0 ] 1 , uS* [ [ u [ 1 , d i m(ul [ 
2111 ]>estretch(u> 


ity(dimCaM2]))''pi,p,l,dim(aM2]iiedim 
_zers(a, zl 

: 


: i i bloq_ Jordan 

: i i din_ Jordan 

C s d i f f s < d i f f s < stretch ( y > > > +d i n_j ord < y > 


C : d i n_ j ord < d i m_zers < a , z "> > +b 1 _.i ord ( a , z > 

MAIN F;fi[i EKACT SEC! 


MAIN F;h[i EKACT SEC: 


: ran k t 0m 5 
: Func 
: Local 0k 
: ref (0m)40m 

:For 0k, rowDim(0ml, 1 , _ 1 

: If sir ing(norm':0mt0l<: 11+0. l/string(0. 1 

: Exit 

: EndFor 

:0k 

: EndFunc 


| t f^lcontro 1 1 1 /Q |Uar |F i nd... |Mode | 

: ua 1 p ( 0n 5 
: Func 

: Local 0p,0z 
: detC0m-0xl-*0p 

: cZeros ( expr ( sir i ng C 0p 1 1 , 0x 1 -+0z 
: EndFunc 


iri-?m)Y F £T Y F3 Y FH t T FE Y Fb^ Y "| 

| ▼ H- |h 1 gebr a |C a 1 c |u t her | Prgn 1 0 |C 1 ean Up | 

■ d i m_ j ord£ d i n_zers( a , z)) + b 1 

_jord(a , zl 



□one 

■ eigUl(m31 

E r r o r s N o n - r e a 1 res u 1 1 

■ eigUl(m31 



4-346e "5 ■ i 3.99996716669 

3. + 7. 6232 1^ 

■ oalp(m31 


f 3 4 J 

■ bl_jord(m3 , 41 

[1 0 

1 0 0 0 0] 

■ b 1 _ j or d( m3 , 31 

[1 1 

0 0 0 0 G] 

bl_.iord<m3 ,3>l 

MAIN RAD EKACT 

SEC! 

13/30 
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Visualizing a Special Line in the 3-D Space 

Nurit Zehavi & David Ben-Chaim 
Weizmann Institute of Science, Israel 


The motivation 

The second author shared with the first author a paper that he had submitted for publication, 
entitled, A Special Point in the 2-D Plane and a Special Line in the 3-D Space. The paper ex- 
plores systems of linear equations that are written in the standard form with coefficients oc- 
curring as consecutive terms of arithmetic sequences. He challenged her: Can we prepare an 
activity for teachers based on this paper using Derive? 

The result 

It is easy to verify, algebraically and graphically, that the ordered pair (-1, 2) satisfies any 
linear equation in two variables that when written in the standard form the coefficients are 
consecutive numbers from an arithmetic sequence.^ 

The solution of the general 3-by-3 system of linear equations with arithmetic sequences for 
coefficients , 

al.x + (al + dl)y + (al + 2dl)z = al + 3dl 

a2.x + (a2 + d2)y + (a2 + 2d2)z = a2 + 3d2 

a3.x + (a3 + d3)y + (a3 + 2d3)z = a3 + 3d3 

is parametric [@1, -2@1 - 1, @1 +2], 

The graphical representation of a linear equation in three variables is a plane. The interpreta- 
tion of the algebraic solution above is that all the planes that represent equations with "arith- 
metic coefficients" intersect in a Special Line (t, -2t — 1, t + 2) (choosing t to be the parame- 
ter). How can we demonstrate this finding applying Derive? 

Derive can place only one surface at a time in the 3D-plot window. However, we superim- 
posed several planes in Acrospin and watched for the intersection line while running Acro- 
spin. The special line in this representation was not very convincing. Coincidentally we re- 
ceived the card from Soft Warehouse with an invitation to "surf the net to their web page. In 
the Derive-Questions-and-Answers section it was suggested to plot the minimum of two sur- 
faces in order to see their intersection. Aha! 

We applied the method to an example taken from Ben-Chaim's paper: 

x + y + z = 1 
x + 7y + 13z = 19 
5x + y - 3z = -7 

The equations were modified to the explicit form z = u(x, y) and the minimum of the three 
planes was plotted (see figure 1). The following figures were plotted by teachers who ex- 
plored the situation. In figure 2 the first equation was replaced by another one with arithmetic 
coefficients, x - z = -2, and the special line appeared again. But do we see the same picture? 
Some teachers studied the 3D plots and compared coefficients, others ran the Acrospin file (it 
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is helpful to plot each of the three planes in a different colour), and a few issued the dot prod- 
uct of the vectors orthogonal to each of the planes and compared the angles between the 
planes. In figure 3, the first equation was replaced by z = y so that the three planes intersect 
in one point (-1, 1, 1). Now a question was raised: Is it possible to see the special line even if 
one of the equations does not fulfil the arithmetic coefficients condition? 



Figure 1 MIN (-x-y+1, - (x+7y-19) ) 


Figure 2 
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The equations , ax + by + cz = d, to consider must be represented by planes that are parallel to 
the identified special line [0, -1, 2] + t[l, -2, 1] (written in vector notation). Thus the dot prod- 
uct [1, -2, l].[a, b, c] = a -2b + c must be zero. In figure 4 such equation, z = x, was chosen to 
replace the middle equation. Not surprisingly, the plot in figure 5 is the same as in figure 1, 
but the coefficients of x + y + z = 10 are not consecutive numbers of an arithmetic sequence. 
In fact, the plane that represents x + y + z= 10 is parallel to one of the planes seen in figure 5. 
In conclusion, the special line activity provides opportunities for teachers to revise and extend 
their knowledge by using the possibilities offered by the software. 



Figure 4 Figure 5 


Nurit and David mentioned in their paper the ACROSPIN animation of their objects. As some of you 
might not know how to combine several DERIVE 3D-plots in different layers and colours in one 
ACROSPIN-file I’ll explain how to do: 

Edit the first plane: 

pll :=-x-y+l. 

Plot the grid in an appropriate colour (Grids 20) in the 3D-plot window and then press Transfer Acro- 
spin Save: MYNAME. Now you may Run Acrospin to view your result so far. Switch back to the Alge- 
bra window, highlight or edit the second plane: 


p!2 :=- (x+7*y-19) /13 
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Plot it again in an appropriate colour. Then Transfer Acrospin Save, confirm MYNAME.ACD pressing 
■ . Skip "Overwrite existing file" pressing N and confirm "Append to existing file" pressing Y . Then go 
back once more to the Algebra window and repeat the same procedure once more with the third 
plane. 


pl3 : = ( 5*x+y+7 ) /3 

Now run Acrospin. 



The right Acrospin animation containing the planes together with the intersection line was produced 
using myACD.EXE tool. (Next DNLM). Josef 


11 E. Sawada, Extending Algebra Concepts with Technology, DNL#12, December 1993 
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Some Notes on DERIVE 2.6 's 
Functions and Limits 

G P Speck, Wanganui, New Zealand 

Using DERIVE 2.6, we have the following results for specific functions, as shown on the printout be- 
longing to these examples: 

(If there are any differences using other - more recent - DERIVE versions I’ve added them. 
Maybe that Al Rich or David Stoutemyer will have some comments on Mr Speck’s paper?) 


Example 1 : F 1 (x) = x SIN(l/x). For this function F 1 (x), all of the DERIVE 2.6 functional values and 
limits agree with the classical results; e.g., 

lim Fl(x) as x — > 0+ = lim Fl(x) as x — > 0- = lim Fl(x) as x — > 0. 

#1: Notati'onDi'gi' ts 6 

#2: Fl(x) := x-SIN^— 

#3: aal := [ f. , 1 - 1 , -1, -0.5, 0, 0.5, 1] 

#4: bbl := VECTOR (F1 (li) , u, aal) 

#5: [ccl := APPEND([x], aal), ddl := APPEND([Fl(x) ] , bbl)] 

#6: [ccl, ddl] 

x L 1 - L -1 -0.5 0 0.5 1 

#7: 

. Fl(x) 1.1752 0.997916 - 0.0833085 - L 0.84147 0.454648 5 0.454648 0.84147 

#8: 1 !mf (! , d) := APPENDCh'm F, STRING O' ) , (x) as x — > 0, d, =) 

#9: eel := [li'mf(l, +) , li'mf(l, -) , li'mf(l, )] 

Edit lim(F1(x),x,0,1) and Edit lim(F1(x),x, 0,-1) 

ffl := rii'm Fl(x) , lim Fl(x) , lim Fl(x)] 

#10: Lx-.0+- X-.0- x-»0 J 

lim Fl(x) as x — > Pi- = 5 

#11: [eel, ffl] 1 = lira Fl(x) as x — > 0- = 5 

. "Hm Fl(x) as x — > 0=5 


Example 2: F2(x) = IF(x = 0, 5, x SlN(l/x)). For this function F2(x), we have the following DERIVE 
2.6 limit results. We see that each of these results DISAGREE with the corresponding classical calcu- 
lus limit results, each of which is 0 if we work with approximate. 

The function F(x):=IF(x=0,0,x SIN( 1/x)) has all limits equal to 0 as x tends to 0, as it is the case for the 
classical limits. 


F2(x) := 

If x = 0 


#12: 

5 

x- SIN (1/x) 


#13: 

bb2 := VECT0R(F2(u) , u 

, aal) 

#14: 

dd2 := APPEND([F2(x) ] , 

bb2) 

#15: 

[ccl, dd2] 


#16: 

f * ‘ 



L F2(x) 1.17520 0.997916 - O.CB33CB5-i. 0.841470 0.454648 5 0.454648 0.841470 












p 14 


G P Speck: Some Notes on DERIVE 's Limits 


D-N-L#23 


#17: ee2 := ["hmf(2, +) , "hmf(2, , "hmf(2, )] 

ff2 := pirn F2(x), Tim F2(x), lim F2(x)"| 

#18: Lx->0+- x->0- x-.0 J 

#19 : [ee2, ff2]’ 

Simplify and Approximate in Version 2.60: 

lim F2(x) as x — > Oh = 5 

#20: lira F2(x) as x — > 0- = 5 

. Inm F2(x) as x — > 0=5. 

Simplify (left) and Approximate (right) in Version 3.10: 




lira F2(x) as x — > Ot = 5 

#23: lira F2(x) as x — > O- = 5 

. lira F2(x) as x — > 0=5 


Example 3 : 


F3(x) = IF(x > 1, I V(x = 0, 0, -1)). For this function F3(x), we have the following DERIVE 2.6 limit 
results. The left and right hand DERIVE 2.60 limits agree with the classical calculus results, but the 
two-sided DERIVE 2.60 limit DOES NOT AGREE with the classical result. 


F3(x) := 

If x > 0 
1 

#24 : If x = 0 

0 

-1 

#25: aa3 := [- 42 , -1, -0.5, 0, 0.5, 1, ^2] 

#26: bb3 := VECTOR (F3(u) , u, aa3) 

#27: [cc3 := APPEND([x], aa3) , dd3 := APPEND([F3(x) ] , bb3)] 

#28: [cc3 , dd3] 

1 1 

x - 42 - 1 . 0 1 42 

#29 : 22 

. F3(x) -1-1-1 0 111 

#30: ee3 [ltmf(3, +) , ltmf(3, -) , ltmf(3, )] 

ff3 ri fra F3(x), lira F3(x), lira F3(x)l 

#31: Lx-.CN- x-.CN X-.0 J 


#32: [ee3, ff3]’ 
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Simplify and Approximate in Versions 2.58 - 3.1 0 and in 6.10: 

"hm F3(x) as x — > Oi- = 0 

#33: Tim F3(x) as x — > 0- = 0 

. lim F3(x) as x — > 0=0. 

Simplify in Versions 3.14: 

lim F3(x) as x — > Oi- = 1 

#34: Tim F3(x) as x — > 0- = -1 

. Tim F3(x) as x — > 0 = ? . 


Example 4: F4(x) = IF(x S 1 N( tt/.t) = 0,1,2). For this function F4(x), the DERIVE 2.6 limits with the 
value 1 DO NOT AGREE with the corresponding classical calculus limits, while those with the value 
2 DO AGREE with the classical limits; also the lim F4(x) as x — > 0 as given below (in approximate 
Mode and if interpreted appropriately) DOES AGREE with the classical calculus limit. 


F4(x) := 

If x-SINCtt/x) = 0 
#35 : 1 

2 

f 111 1 1 12221 


#36: aa4 := -2, 


#37: bb4 := VECTOR (F4(u) , u, aa4) 


1000000 1000000 10 3 9 27 


#38: [cc4 APPEND([x] , aa4) , dd4 := APPEND([F4(x) ] , bb4)] 

#39: [cc4, dd4] 

no problems with DERIVE 6.10: 


12 2 


1000000 1000000 10 3 9 27 


F4(x) 2 2 2 


2 2 2 


-6 -6 

x -i -2 -0.5 -0.333333 -0.25 - 10 0 10 0.1 0.666666 0.222222 0. 


0740740 


L F4(x) 2 2 1 1 1 

there were problems in 1 996 with DERIVE 3.14: 


1111 


This is from 1996 


-i -2 - 


1000000 


2 2 


1 0-SIN - 


1000000 10 


>, 1, 2 










#42: F5(x) := |SIGN(x)| 

#43: aaB := [ i. , 1 - i, -1, -0.5, 0, 0.5, 1] 

#44: bbB := VECT0R(F5(u) , u, aa5) 

#45: [cc5 := APPEND([x], aa5) , dd5 := APPEND( [F5 (x) ] , bb5)] 

#46: [cc5, dd5] 

x lI-l -1-0.5 00.51' 

#47: 

. F5(x) 11 11111. 

#48: eeB := [1n'mf(5, +) , In'mfCB, -) , T n rm f ( 5 , )] 

ffB := Him F5(x), Inm F5(x), Inn F5(x)l 

#49 : Lx -Of x-0- x-0 J 

Tim F5(x) as x — > Of = 1 

#50: [eeB, ffS] 1 = Tim F5(x) as x — > 0- = 1 


Tim F5(x) as x — > 0 


1 
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The Derivative: 

#51: [F5' (x) , F5 r (0)] = [0, 0] 

This is from the Utility file NurnericalApproximation.rnth: 

DIF_NUMERIC (y , x, xO, h, n) := 

If NUMBER (n) 

#52 E / (-1) A k ■ COM B / n , k)-LIM/y, x, xO - n-h/2+k-h), k, 0, n)//-h) A n 

DIF_NUMERlC(y , x, xO, h, 1) 

DIF_NUMERIC/y , x, xO, h, 1) 

#53: DIF_NUMERIC(F5(x) , x, 0, h) = 0 


Example 6: 

F6(x, v) = IF(y = x 2 , 1 , 0). Each of the first two DERIVE 2.60 limits below DISAGREES with the cor- 
responding classical limit for a function of two variables, since the classical calculus limit fails to exist 
in each case. The third limit below AGREES with the corresponding classical calculus one. 

The same co mm ents can made be here for the LIM2 results as were made for the LIM results above. 
Note that neither the DERIVE 2.60 LIM nor LIM2 (being limits along rays through the point in ques- 
tion) are designed to handle the classical limit of a function of the F6 type, even in the absence of the 
difficulties shown here. 



#56: APPEND/ [[n, 

F/cl) , F/al) , F/bl)]], 

VECTOR /[n , F/cl), F/al), 

F/txL)], n, 

2, 7)) 

APPEND/ [[x, 

y, F6/X,y)]], VECTOR/r 

aa6 , aa6 , F6/aa6 , 

aa6 ) ] , 

k, DIM/aa6))) 

#57: 

l 

k ,1 k , 2 k ,1 

k ,2 J 



x y F6/x,y) 



i 

-1 

i 


L 

L 

0 


1 

1 

1 

#58: 

0 

0 

1 


1 




— 

1 

0 


2 




. -1 

-2 

0 

ee6a := 

[Tim 

F6/x,y) 

#59: 


L 


ee6b := 

film 

F6/x,y) 

#60: 


L 


ee6c := 

Mini 

F6/x,y) 


as 


as 


as 


[x,y] — > 
[x,y] — > 
[x.y] — > 


[0,0] = 

[x. 

Tim 

y]-*[o, 

0] 

F6/x, 

y) 

[1,1] = 

[x, 

Tim 

y ]->[l, 

i] 

F6/x , 

y) 


[ 1 / 2 , 1 / 2 ] = , 


# 61 : 


lira F6(x, y) 

[x, y Ml/2, 1/2] 
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Version 3.14 

Tim F6(x,y) as [x,y] — > [0,0] = 0 

#63: lim F6(x,y) as [x,y] — > [1,1] = 0 

. Tim F6(x,y) as [x,y] — > [1/2 ,1/2] = 0 

Version 2.6 and 6.10 

Tim F6(x,y) as [x,y] — > [0,0] = 1 

#64: lim F6(x,y) as [x,y] — > [1,1] = 1 

. "h'm F6(x,y) as [x,y] — > [1/2 ,1/2] = 0 

LIM2(u, x, y, xO, yO) := lira 


#65: 




[y 

, X] 

-•[yO + @1- (x ■ 

#66: 

ee6d := 

[’ll m2 F6(x,y) 

as 

[x 

,y] 

--> [0,0] = , 

#67: 

ee6e := 

[11(112 F6(x,y) 

as 

[x 

,y] 

-> [1,1] = , 

#68: 

ee6f := 

1im2 F6(x,y) 

as 

[x 

,y] 

-> [1/2, 1/2] 

#69: 

[ee6d , 

ee6e, ee6f] 





Version 3.14 







’ll m2 F6(x,y) as 

[x,I 

/] 

— > 

[0,0] = 0 

#70: 

"h'm 2 F6(x,y) as 

[x ,: 

/] • 

— > 

[1,1] = 0 


L 1 l m2 F6(x,y) as [x,y] — > [1/2 ,1/2] = 0 

Version 2.6 and 6.10 



"Hm2 F6(x,y) as [x,y] — > [0,0] = 1 

#71: In m2 F6(x,y) as [x,y] — > [1,1] = 1 

. 1 n m2 F6(x,y) as [x,y] — > [1/2, 1/2] = 0 


A number of the difficulties illustrated in the preceding examples can be traced to the way in which 
DERIVE 2.60 limits behave relative to functions using the IF( ) statement - however, as Example 5 
shows, this is not the only difficulty. The way in which DERIVE 2.60 limits react relative to the IF( ) 
statement functions given above can be determined without too much analysis; but other more compli- 
cated multiple IF( ) statement functions can be displayed where it is not nearly so easy to predict just 
what the DERIVE 2.60 limit result will be. 

I have created a variety of simple DERIVE "programs" to bypass some of the one and two variable 
limit difficulties, but have not come up with a totally satisfactory solution to the inherent limit prob- 
lems. Example 5 above illustrates a particularly annoying problem in that DERIVE returns a numerical 
value for the composition of two pre-defmed functions at an independent variable value which 
"shouldn’t" be in the domain of the composition. Flow many other examples of this kind are preva- 
lent? None? Several? 

If we were confident that DERIVE will produce appropriate functional values for all predefined and 
author defined (using IF( ) where needed) functions and compositions of such functions, then we 
would be in a position to seriously tackle the resolution of the DERIVE limit difficulties. Notwith- 
standing the difficulties raised herein, for the very great majority of DERIVE functions of whatever 
type that I have used, correct functional values are returned. As long as this is the case, DERIVE does 
provide two very powerful tools for investigating important features, including the limit, of non-trivial 
functions: the easy display of a table of functional values and functional plots. As an illustration, Ex- 
ample 7 in the printout that follows shows the verification (by way of a table) of an important feature 
of the very interesting non-analytic function 
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Example 7 : An example of a non- elementary type of local minimum 


F(x) = A(x) + 5(x)sin 

C(x) + C f- 
A(X) = - W 


where 


C(x) - C 


and B(x) = 


C(x) = e * 2 for x * 0 and 

F(x) = 0 for x = 0. 


C(x) + C — C(x) - C — 

#72: L 2 J l 2 

A(x) := , B(x) := 

2 2 


#73: - 1/x 

C(x) s 


#74: F(x) := A(x) + B(x)-SIN — 


f r 1 1 f ( 1 

2 SIN 2 SIN — 

#7S : - 1/x L x J 1 - 4/x 1 Lx 

F(x) := e ■ + + e ■ - 

L 2 2 J {. 2 2 


+ (2- n - 1) -7T 

2 


+ (2- n + 1) -n 

2 


■0.0005 

■0.0004 

■0.0003 

■ 0.0002 

■ 0.0001 


1-0. fi-0. 7-0. 6-0. 5-0. 4-0. 3-0. 2-0.1 0.1 0.2 0.3 0.4 0.5 0.6 

1 - 0.0001 


#77: [PrecisionDigits 20, Notation Dig its := S ] 

#78: APPEND([[n , al, hi, cl]], VECT0R([n, al, hi, cl], n, 2, 7)) 


n 

al 

bl 


cl 

2 

0.0909456 

0 , 0707355 

0. 

0578745 

3 

0. 0378745 

0.0489707 

0. 

0424413 

4 

0.0424413 

0.0374482 

0. 

0335063 

5 

0. 0335063 

0.0303152 

0. 

0276791 

6 

0.0276791 

0.0254647 

0. 

0235785 

7 

0. 0235785 

0.0219524 

0. 

0205361 
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#80: APPEND [[n, SIN(l/al), SIN(l/bl), SIN(l/cl)]], VECTOR n 



' n SINCl/al) SIN(l/bl) SIN (1/ cl) ' 

2-1 1 -1 

3-1 1 -1 

#81: 4-1 1 -1 

5- 1 1 -1 

6 - 1 1 -1 

.7-1 1 -1 

#82: APPEND([[n , F(cl) , F(al), F(tl)]], VECT0R([n, F(cl) , F(al), F(tl)], n, 2, 7)) 


n 

F(cl) 

F(al) 

F(ti) 


-S19 

-211 

-87 

2 

2.26943-10 

9.34497-10 

1.59267-10 


-96 S 

-519 

-182 

3 

3.79477-10 

2.26943-10 

8.00207-10 


-1S48 

-965 

-310 

4 

4.369-10 

3.79477-10 

2.05949-10 


-2268 

-1548 

-473 

5 

3.46344-10 

4.369-10 

2.71519-10 


-3125 

-2268 

-670 

6 

1.89043-10 

3.46344-10 

1.83367-10 


-4120 

-3125 

-902 

7 

7.10467-10 

1.89043-10 

6.34348-10 



It can be proved for our F(x) defined above for all real x: 

F(x) has derivatives of all orders at each real x; 

F(x) has a relative minimum at x = 0; 

al > bl > cl and F(cl) < F(al) < F(bl) for n > 2. 

Thus, F'(x) is not > 0 throughout any interval (0,b), nor is F'(x) < 0 throughout any interval (a,0). 

The DERIVE matrices displayed above provide simple verification that al > bl >cl and F(cl) < F(al) 
< F(bl) for 2 < n < 7. Such simple DERIVE verification is very handy indeed in providing confidence 
in one’s theoretical results for "intricate" functions such as our F(x) given here. 

The following is a "graphical" view of the relative magnitudes of al, bl, cl and F(al), F(b 1), F(cl). 


F 


for n >= 2: 


cl < bl < al 
F(cl) < F(al) < F(bl) 
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A word of caution! If we define F(x) in the following manner without Simplifying F(x) as was done 
above, then the matrix of values for F(al), F(bl) and F(cl) may be totally inaccurate!! 


#98 : 


F(x) := A(x) + B(x) ■ SIN 



Approximated value table without simplifying F(x): 

#99: APPEND([[n, F(cl), F(al), F(tl)]], VECT0R([n, F(cl) , F(al), F(bl)], 


n 

F(cl) 

F(al) 

F(ti) 


-260 

-106 

-87 

2 

2.38192-10 

4.83347-10 

1.59267-10 


-483 

-260 

-182 

3 

3.08 008-10 

2.38192-10 

8.00207-10 


-774 

-483 

-310 

#100: 4 

1.0451-10 

3.08008-10 

2.05949-10 



-774 

-473 

5 

0 

1. 0451-10 

2.71519-10 




-670 

6 

0 

0 

1.83367-10 




-902 

. 7 

0 

0 

6.34348-10 


2, 7)) 


We see that the F(cl) and F(al) values in this matrix are totally different from the correct function 
values in the simplified F(x) matrix above, while the F(bl) values agree for the two matrices. A classi- 
cal problem in numerical analysis underlies the difficulty here: 

For our function F(x) e A (-4/x A 2) is so much smaller than e A (-l/x A 2) for al with n = 2, for example, 
that the evaluation of A(al) with PrecisionDigits: 20 will give the value A(al) = (4.83347 10 A -106)/2 
and similarly for B(al). If we used PrecisionDigits: 100 in this case the following shows that the 
DERIVE F(al) value would be the correct value: 


#101: n := 2 

#102: [A(al), B(al)] 

-53 -53' 

#103: Ll. 55458-10 , 1.55458-10 

#104: F (al) 

-211 


#105: 9.34497-10 
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The Trigonometric 
Super Box 


Josef Bohm, Wiirmla, Austria 


In DNL#9 from March 1993 I presented a "Black Box" to solve oblique triangles. At the occasion 
preparing a workshop for the DERIVE Conference 1 996 I tried to extent this "Black Box" to a real 
"Superbox" in order to build a "Circuit of Boxes" for this special subject. 



For all functions and formulae we define a 
triangle with common names for sides and 
angles: 

(I based my workshop on the problem: 

AABC(BC = 120, AC = 80, ZABC = 20°)). 

This is the most difficult case, because we 
have to consider two solutions! ! 


1 The White Box 

In the White Box phase we would derive the Sine Rule and the Cosine Rule. The scientific pocket 
calculator is sufficient support at that time. The first step using DERIVE could be editing the Cosine 
Rule and the Sine Rule. The students must know which of the both rules applies for solving a problem, 
which variables can be substituted and then let DERIVE solve the respective equation. I’d like to call 
this position as 

2 Not a Box at all 

The start of solving the given problem could look like as follows: 

Not a Box at All 

User : 

#1: InputMode := Word 

User : 

2 2 2 

#2: s3 = si + s2 - 2 • si • s2 • C0S(a3) 

User : 

#3: sl-SIN(a2) = s2-SIN(al) 

User : 

#4: 80 • SIN(ot°) = 120-SIN(20°) 

Sol ve(#4,a) : 

#5: SOLVE(80 • SIN(oi°) = 120 • SIN(20°) , <x, Real) 

Approx(#5) : 

#6: a = -210.8658824 v a = 149.1341175 v « = 30.86588247 
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Take care that you set Manage Trigonometry Angle: Radian (in DERIVE 6: Options > Mode Settings 
> Angular unit > Radians). So you have to substitute for the angle a2 in line #3 a deg. In the next 
phase we provide four functions which give the solutions of the equations in lines #2 and #3 for the 
angle or the side. This package of functions represents the 


3 The "Gray" Box 


cs_a3(side_l, 

side_2, angle_3) 

returns side_3 

ca_sl(si de_l , 

side_2, side_3) 

returns angle 1 

ss_a2(si de_l , 

angle_l, angle_2) 

returns side_2 

sa_s2(si de_l , 

side_2, angle_l, k) 

returns angle 2 (k = 1 => angle 2 < 90°, 



k = 2 => 90° < angle_2 < 180°, 



k = 1 by default.) 


The students must still know which of the rules applies. They also have to consider the case with the 
possibility of obtaining two solutions. To explain the "names" of the functions: 

cs_a3 = use the Cosine Rule to find the Side opposite to Angle 3 . 

(You will find the "Gray Box functions" in the DERIVE printout at the end of the contribution.) 

4 The "Black Box" 

These four functions from the "Gray Box" can be developed together with the students. The next ones 
are more difficult. We use the Gray Box Functions and combine them to a real "Black Box" to solve 
the triangles in one step. It will depend on the students' skill in working with DERIVE to produce the 
next package. This could be much easier if the maths teacher is also the information technology 
teacher. We can build in some control structures (considering the triangle’s inequality, considering 
one, two or none solution in case of SSA, ....). I could imagine to produce the most difficult module 
together with the students and let them produce the remaining ones in groups by their own. 

Last year I was short of time teaching trigonometry. I only had time to deal with the trigonometric 
functions and their applications in connection with right angled triangles. But I didn’t want to miss the 
applications of solving oblique triangles in plane and space. I wanted the students to use strategies for 
calculating polygons and problems in 3D-space. So I provided this Black Box in a Utility file and an 
instruction similar to the following: 

Load the Utility file TRIGO.mth. You then will have 
four functions available for solving triangles. The 
provided functions are: 

sss(side_l, side_2, side_3) 
sas(side_l, angle_3, side_2) 
asa(angle_2, side_l, angle_3) 
ssa(side_l, side_2, angle_l) 


(/ don ’t want to repeat that part from DNL#9. 1 use these functions for building the SUPER BOX.) 
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But according to Bernhard Kutzler [2,3] the real Super Box would be a function like 

triangle ( side 1, side_2, side_3, angle 1, angle 2, angle 3) 
with three parameters known and three u nk nown. DERIVE should do all the remaining work: check if 
the triangle really exists, decide the case and calculate the missing sides and angles. Then our example 
should look like that: 

triangle ( 12 0 , 80, side_ab, angle_bac, 20°, angle_abc) . 

I faced the challenge to develop such a function and we reproduced the ideas in the workshop. 


5 The "Super Box" 

Depending on the u nk nown parameters in TRIANGLE (...) DERIVE has to decide which of the Black 
Box functions applies. Additionally we want to implement the until now untreated case of SAA. (This 
is new even for the participants of the Workshop in Bonn). In the Black Box phase we would calculate 
the third angle by hand and then apply ASA. That is now the moment when not only the - not too dif- 
ficult - implementing of mathematical rules will take place but when you also - perhaps with your 
students - have to consider an algorithm. In the following I will present one possibility (I didn’t find a 
better one): 

I set up a table with 6 columns for the parameters and sufficient rows to represent all the possible 


combinations to define a triangle from SSS to AAA - (next page). As j^J = 20 we need 20 rows. Sur- 
prisingly a pattern emerges. This pattern is the key to solve the problem. The parameters have to be 
edited as a vector v and the existence of values for its elements creates the code for the application of 
well known tools. While C2 (v) picks out the SAA-case, CODE (v) gives a unique relation between 
the code numbers and the remaining problems (SSS, ASA, AAA and SSA). The IF-construct returns 
the value 0 if DERIVE cannot decide whether the expression is 0 or not. And that is the case, if the 
expression is a not predesignated variable or a given constant value. 

c2 (v) : =IF (v SUB l*v SUB 4 *v SUB 5>=0 OR (v SUB l*v SUB 4*v SUB 6>=0 OR 
(v SUB 2 *v SUB 4 *v SUB 5>=0 OR (v SUB 2*v SUB 5*v SUB 6>=0 OR 
(v SUB 3*v SUB 4 *v SUB 6>=0 OR v SUB 3*v SUB 5*v SUB 6>=0) ) ) ) , 1 , 0, 0) 

code (v) :=SUM (VECTOR (IF (v SUB i=0, i, i, 0) , i, 6) ) 

For example see some codes: 

code([102, 87, 113, a, [3 , *]) = 6 
code([x , y, z, 45°, 65°, 70°]) = 15 

code([120, 80, 5ide_AB, ang1e_BAC, 20°, ang1e_ABC]) = 8 


For the cases SAS and ASA only a permutation of the parameters is necessary. AAA is special: 
DERIVE returns the proportions of the sides (see one example!). The most difficult problem occurs in 
the cases of SSA and SAA, because there is no unique relation between the code numbers and the 
given sides and angles. But as you can see in the special printed numbers in the table there is an 
unique incidence of pairs of numbers. We can use the fact that the product of these pairs is defined and 
greater 0 to identify the problem. See first the table followed by some examples : 
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side_1 = v-i 

side_2 = v 2 

side_3 = v 3 

angle_1=a 1 

angle_2=a 2 

angle_3=a 3 

Z = Code 

Case 

© 



4 

© 


10 

SAA 

© 



4 


© 

11 

SAA 


© 


© 

5 


11 

SAA 


© 



5 

© 

13 

SAA 



© 

© 


6 

13 

SAA 



© 


© 

6 

14 

SAA 

1 

2 

3 




6 

SSS 

1 

2 




6 

9 

SAS 

1 


3 

5 



9 

SAS 


2 

3 

4 



9 

SAS 

1 




5 

6 

12 

ASA 


2 


4 


6 

12 

ASA 



3 

4 

5 


12 

ASA 




4 

5 

6 

15 

AAA 

1 

© 


© 



7 

SSA 

© 

2 



© 


8 

SSA 


2 

© 


© 


10 

SSA 


© 

3 



© 

11 

SSA 

1 


© 

© 



8 

SSA 

© 


3 



© 

10 

SSA 


#56: 


#57: 

#58: 


#59: 


triang1e([side_a, side_b, 102, 118.34°, 8, 41.26°]) 
IT s i de_b = 53.9129 ' 
side_a = 136.130 
. 8 = 20.4 degrees . 

triangle([80, 120, PR, a_QRP, 20°, a_PQR]) 
so!ution_number = 1 
a_QRP = 13.1801 degrees 
a_PQR = 146.819 degrees 


PR = 192.014 


#60: triang1e([120 , 80, PR, a_QRP, 46°, a_PQR]) = no_triang1e 
#61: triangle([120 , 80, PR, a_QRP, 20°, a_PQR]) 



so1ution_number = 1 


solution_number = 2 



a_QRP = 30.8658 degrees 


a^_QRP = 149.134 degrees 



a_PQR = 129.134 degrees 

» 

a^PQR = 10.8658 degrees 



PR = 181.432 


PR = 44.0934 



# 62 : 
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#63: triangle([a, b, c, 30°, 60°, 90°]) 


#64: simplified 


#65: 


only similar triangles are possible: 
a = tAb=-/3t*c=2t 


#66: approximated 


#67: 


only similar triangles are possible: 
a = t a b = 1.73205 t a c = 2 t 



CIRCUIT OF BOXES 

The workshop and, of course, this contribution 
should demonstrate that you are able to "pro- 
gram" even in DERIVE. You can use control 
structures and "error catching". 

There are many possibilities to change from one 
box phase to another one. I try to illustrate that 
fact. Y ou can step into this circuit at any position 
and then move to all directions. 


Now, in 2010 with DERIVE 6.10 and its programming features available I would doubtless try 
writing a compact DERIVE program. There are some differences to the code given in 1996 
because of changes in DERIVE's functionality (e.g. SOLVE in 1996 and 
SOLVE/SOLUTIONS in later versions). In later versions we can differ between ASIN and 
ARCSIN. This is also not considered in this package. But the package works as it is. 

#1: SUPERBOX. MTH - The Trigonometric Super Box 

#2: [InputMode := Word, CaseMode := Sensitive, Angle := Radian] 

#3: [TimesOperator Implicit, Notation :r Decimal, NotationDigits 6] 

#4: c_angle(al, a2) 7r - al - a2 

tri_inequ(sl, s2, s3) := 

If si + s2 < s3 v (s2 + s3 < si v si + s3 < s2) 

#5: 1 

0 


#6: 


The Gray Box Functions: 


2 2 

#7: cs_a3(sl, s2, a3) := ./(si + s2 - 2 si s2 C0S(a3)) 


ACOS 


# 8 : 


f 2 2 2 

s2 + s3 - si 

2 s2 s3 


ca_sl(sl, s2, s3) 


si SIN(a2) 

#9: ss_a2(sl, al, a2) := 

SIN(al) 


# 10 : 


ASIN 


s2 SIN(al) 


si 


(3 - 2 k) + Jr (k - 1) 


sa_s2(sl, s2, al, k := 1) := 


1 ' 
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|The Black Box| 


#12: We know the three sides si, s2, s3. 

sssCsl. s2, s3) := 

If tri inequ(sl, s2, s3) = 1 
#13: no_tri angle 

[angle_opposite_sl = si, ca_sl(sl, s2, s3) degrees; angle_opposite_s2 = s2, ca_sl(s2, s3, si) degrees; 


#14: 


#15: 


#16: 


#17: 


#18 

#19 

#20 


# 21 : 


# 22 : 


#23: 


#24: 

#25: 

#26: 


#27: 


#28: 


angle_opposite_s3 = s3, ca_sl(s3, s2, si) degrees] 

We know two sides si and s2 and the angle a3 formed by them. 

a3 

side_opposite_a3 = degrees : cs_a3(sl, s2, a3) 

1 ° 


sas(sl, a3, s2) := 


angle_opposite_sl = si 
angle_opposite_s2 = s2 

We know side s3 and both adjacent angles al and a2. 


ca_sl(sl, s2, cs_a3(sl, s2, a3)) degrees 
ca_sl(s2, si, cs_a3(sl, s2, a3)) degrees 


angle opposite. s3 = s3 


asaCal, s3. a2) ;= 


al 

side opposite al = degrees 

1 ° 


c_angle(al. a2) 
1 ° 


degrees 


ss a2(s3, c angle(al, a2). al) 


a2 

side opposite a2 = degrees ss a2(s3, c angle(al, a2), a2) 

1 ° 

sides si and s2 and an opposite angle are given 

One solution, no solution or two solutions are possible! 

I need an auxiliary function: 


sol ution_of_tri angle = k 


angle_opposite_s2 = s2 


sa_s2(sl, s2, al, k) degrees 


ssah(sl, s2, al, k) := 


thi rd_angle 


c_angle(al, sa_s2(sl, s2, al, k) 1°) 
l e 


degrees 


third_side : ss_a2(sl, al, c_angle(al, sa_s2(sl, s2, al, k) 1°)) 

Decision of cases using the auxiliary function 

ssa(sl, s2, al) := 

If si > s2 

ssah(sl, s2, al, 1) 

If s2 SIN(al)/sl > 1 
no triangle 

VECTOR([ssah(sl, s2, al, k)], k, 1, 2) 

Some Examples: 
sss(102, 61, 109) 

angle opposite si = 102 : 66.9903 degrees 

angle opposite s2 = 61 : 33.3984 degrees 

. angle_opposite_s3 = 109 : 79.6111 degrees . 

solution of triangle = 1 


ssa(100, 160, 33 1°) 


angle_opposite_s2 = 160 
thi rd_angle 
thirtLside 


60.6241 degrees 
86.3758 degrees 
183.24 


sol ution_of_tri angle = 2 


sas(85 .3. 58.6 1°, 23.4) = 


angle opposite s2 = 160 
third angle 
third_side 

side_opposite_a3 =58.6 degrees 
angle_opposite_sl =85.3 
angle opposite_s2 = 23.4 


119.375 degrees 
27.6241 degrees 
85.1343 


75.7876 

106.119 degrees 
15.2802 degrees 


#29: 


ssa(100, 160, 53 1") = no triangle 
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so1ution_of_triangTe = 1 


#30: 554(160, 100, 33 1”) = 


#31: asa(67.3°, 2987, 50.1°) = 


ang1e_opposite_s2 - 100 
thi rcLangle 
third_side 

ang1e_opposite_s3 = 2987 
side_opposite_al =67.3 degrees 
side_opposite_a2 =50.1 degrees 


19.9012 degrees 
127.098 degrees 
234.312 


62.6 degrees 
3103.82 
2581.07 


#32: 

#33: 


The Super Black Box 


code(v) := I(VECTOR(IF(v =0, i, i, 0). i, 6)) 
#34: i 


#35: 


#36: 


c2(v) := 

If Vil Vi4 v;5 £ 0 v (v;l vi4 vi6 z 0 v (v;2 Vi4 v;5 z 0 v (v;2 v;5 vi6 £ 0 v (vi3 vi4 v*6 £ 0 v vi3 Vi5 Vi6 t 0)))) 
1 
0 
0 

sssx(sl, s2, s3, al, a2, a3) := 

If tri_inequ(sl, s2. s3) = 1 
no_triangle 

[LIM(a_, a_, al) = ca_sl(sl, s2, s3) degrees; LIM(a_, a_, a2) = ca_sl(s2, s3, si) degrees; LIM(a_, a_, a3) = 


ca_sl(s3, s2, si) degrees] 

s_s_s(v) := sssx(v , v , v , v , v , v ) 
#37: 1 2 3 4 5 6 


lim s_ = cs_a3(sl, s2, a3) 
s_-»s3 


#38: sasx(sl, s2, s3, al, a2, a3) := 


lim a_ = ca_sl(sl, s2, cs_a3(sl, s2, a3)) degrees 
a_-*al 

lim a_ = ca_sl(s2, si, cs_a3(sl, s2, a3)) degrees 
a_-a2 


s_a_s(v) := 

If vi4 > 0 

sasx(vi2, vi3, vil, v;5, Vi6, v;4) 

0 

#39: If vi5 > 0 

sasxCvil, vi3. vi2, Vi4. va 6. Vi5) 
0 

sasxCvil, vi2, vi3, Vi4, vi5, vi6) 


#40: asax(sl, s2, s3, al, a2, a3) := 


c_angle(al, a2) 

lim a_ = degrees 

a_-a3 1° 

lim s_ = ss_a2(s3, c_angle(al, a2), al) 
s_-sl 

lim s_ = ss_a2(s3, c_angle(al, a2), a2) 
s_-s2 


a_s_a(v) := 

If vil > 0 

asax(vi2, Vi3, vil. ViS, vt6, vi4) 

0 

#41: If Vi2 > 0 

asaxCvxl, Vi3, vi2. vj4, vi6. Vi5) 
0 

asaxCvil, vi2, vi3, vi4, ViS, vi6) 


#42: ssahx(sl, s2, s3, al, a2, a3, k) 


solution_number = k 

lim a_ = sa_s2(sl, s2, al, k) degrees 
a_-»a2 

c_angle(al, sa_s2(sl, s2, al, k) 1°) 

lim a_ s degrees 

a_-*a3 1° 

lim s_ = ss_a2(sl, al, c_angle(al, sa_s2(sl, s2, al, k) 1°)) 
s_-»s3 


ssax(sl, s2, s3. al, a2. a3) := 

If si * s2 

ssahxCsl, s2. s3, al. a2. a3. 1) 

If s2 5IN(al)/sl > 1 
no_triangle 

(C0PR03ECTI0N(VECT0R([ ssahx(sl , s2, s3, al, a2. a3, k) ] . k, 1, 2 )))a1 


#43: 
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#44: 


aux(v, zl, z2, z3) := ssax(v , v , v , v , v , v ) 

zl z2 z3 zl + 3 z2 + 3 z3 + 3 


s_s_a(v) s 

If Vi2 Vi4 > 0 
aux(v. 1. 2, 3) 

0 

If vil Vi5 > 0 
aux(v, 2, 1, 3) 

0 

If V|3 vi4 a 0 

#45: aux(v, 1, 3, 2) 

0 

If vjl Vi6 > 0 
aux(v, 3. 1, 2) 

0 

If Vi3 vjS i 0 
aux(v, 2, 3, 1) 

0 

aux(v, 3, 2, 1) 

sols(v) := SOLVEtfv = t. v SIN(v ) = v SINCv ). v SINCv ) = v SINCv )1. [v . v . v 1) 
#46: [ 1 1 5 2 42 6 3 5 J L 1 2 3j 


special(v) s 

If v»4 + V|5 + vi6 * 180 1° 

#47: "Sum of angles incorrect!" 

["only similar triangles are possible:"; (sols(v))il] 


lim s_ = ss_a2(sl, al, c_angle(al, a2)) 
s_-<s3 


#48: saax(sl, s2, S3, al, a2, a3) :; 


lim s_ = ss_a2(sl, al, a2) 
s_-<s2 


c_angle(al, a2) 

lim a_ = degrees 

a_->a3 1° 


aux2(v, zl, z2, z3) :; saax(v , v , v , v , v , v ) 

#49: zl z2 z3 zl + 3 z2 + 3 z3 + 3 


s_a_a(v) :r 

If Vil Vi5 > 0 

aux2(v, 1, 2, 3) 

0 

If Vil Vi6 > 0 

aux2(v, 1, 3, 2) 

0 

If Vi2 Vi4 > 0 

#50: aux2(v, 2, 1, 3) 

0 

If Vi2 Vi6 > 0 

aux2(v, 2, 3, 1) 

0 

If Vi3 Vi4 > 0 

aux2(v, 3, 1, 2) 
0 

aux2(v, 3, 2, 1) 

triangle(v) := 

If c2(v) = 1 
s_a_a(v) 

If code(v) = 15 
special (v) 

If code(v) = 6 
#51: s_s_s(v) 

If code(v) = 9 
s_a_s(v) 

If code(v) = 12 
a_s_a(v) 
s_s_a(v) 


#52: 


Examples : 


#53: triangle([side_a, side_b, 102, 118.34°, (5, 41.26°]) 


side_b = 53.9128 
s i de_a = 136.130 


#54: triangle([80, 120, PR, a_QRP, 20°, a_PQR]) 


(5 = 20.4 degrees 
solution_number = 1 
a_QRP = 13.1801 degrees 
a_PQR = 146.819 degrees 


PR = 192.014 


#55: triangle([120, 80, PR, a_QRP, 46°, a_PQR]) = no_triangle 
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#56: triang1e([120, 80, PR, a_QRP, 20°, a_PQR]) 


#57: 

#58: 

#59: 

#60: 

#61: 


solution_number = 1 
a_QRP = 30.8658 degrees 
a_PQR = 129.134 degrees 
PR = 181.432 
triangle([a, b, c, 30°, 60°, 90°]) 
simpl if led 

only similar triangles are possible: 
a : t A b i /3 t A C : 2 t 
approximated 

only similar triangles are possible: 
a = t a b = 1.73205 t a c = 2 t 


solution_number = 2 
a_QRP = 149.134 degrees 
a_PQR = 10.8658 degrees 
PR = 44.0934 


tri angl e(v) is the main function. According to the code numbers it assigns each problem to its solu- 
tion procedure. 

There is one more improvement possible. I will try to append the feature to plot the triangles. I have in 
mind a function PLOT(v) which returns the coordinates of the vertices. One vertex will be in the ori- 
gin and one side will lie on the x-axe. Maybe you can find it on the file on the diskette. 

I’ll close with the same example which was worked through using the Black Box in DNL#9. The prob- 
lem was to solve a Pentagon ABODE. Do it once more testing the Super Box: 

(AB = 30, BC = 52, AE = 30, AD = 65, BD = 60, EC = 68, Z.EAB = 106.1°; ED = ?, CD = ? 

Solution: ED = 45.8171; CD = 36.31 14) 


jr^Tfll gebra|Calc |Other|Prgn 10 |C tear 


a-z... 


2 (3 cost 6] + 4 -since)) * rl(6] Done 

sss(102, 223, 158) 

[24.172773 116.458597 39.368630] 

ssaCSO, 120,20] 

T 30. 865882 129.134118 181.4327531 

L 149. 134118 10.865882 44.093471 J 

sas(23. 56 , 102.69, 18.77] 

[ 43 ■ 826820 33.191525 33 . 4331 80] 


asa<52 -54.17-5 .78 - 6 >1 


trigo\trigo[ 105.33,82.17 

,ab, 11 0.1 3, £y) 

105.33 43.4302 82.17 

110,13 22.7765 47.0935 




w_beta,w_gamma ) 

seite_c 

82.17 

w_beta 

22.7765 

w_gamma 

47.0935J 


[1] The DERIVE News Letter #9, March 1993, pl5-pl7 

[2] Mathematik unterrichten mit DERIVE, Bernhard Kutzler, Addison-Wesley 1995, ISBN 3 89319 860 1 

[3] Improving Mathematics with DERIVE , Bernhard Kutzler, Chartwell-Bratt 1996, ISBN 0 86238 422 2 


As you can see on the screen dump I have tried successfully to implement the "Black Box" 
on the TI-92. I will go further to implement the Superbox, too. Maybe you can find it on the 
diskette of '96. 

Comment in 2010: 1 did not develop the Super Box for the TI-92. But I can offer the SUPER 
BOX for the Tl-Nspire in combination with a short introduction in programming the Tl-Nspire 
in combination with developing a “TRIGO-Library" for the next DNL#77. Josef 
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DREIECK.MTH - 
-TRIANGLE.MTH (2) 

Bernhard Wadsack, Vienna, Austria 

This is the second part of Bernhard's article to plot a triangle together with its remarkable points and circles. Bernhard sent a 
paper which he gives his pupils to work with this Utility file. I add this paper at the end of the contribution. 

You should set the following Options: 

Options Input Word, Plot Options State Connected, Options Auto Change Color No 

#1: TRIANGLE. mth should be preloaded - as utility file 

#2: L0AD(D: \DfD\DNL\DNL96\MTH23\TRIANGLE . mth) 

#3: Die Winkel - the angles 

(2 2 2 1 

ISO lb + 1c - la 

#4: a \- ACOS 

rr L 2 ■ lb ■ 1 c J 

( 2 2 2 1 

ISO la + 1c - lb 

#5: |J := ACOS 

7T C. 2- la- 1c ^ 

( 2 2 2 1 

ISO la +lb -1c 

#6: y := ACOS 

7T ^ 2- la- lb j 

#7: Langen der Hohen - Lengths of the altitudes 

#8: [lh_a := lcSIN(p°), lh_b := la. SIN(y°) , lh_c := lb-SINCO] 

#9: Die Ankreismittelpunkte - Centers of the adj, circles Al_, A2_, A3_ 


A3_ := x3 + 


Cxi - x3) TAN Kyi - y3) 


Cyl - y3) TAN . (x3 - xl) 

2 2 l 2 J 


A2_ := x2 + 


Al_ := xl + 



lb 


lb 


u 


U 1 



— 

■ (x3 - x2) 

TAN 

Ky3 

- y2) 

2 


2 1 

c 2 J 



la 


la 


u 


U 1 



— 

■ (x2 - xl) 

TAN 

Ky2 

- yD 

2 


2 1 

c 2 J 




lb 


lb 


U 


U 1 

r 'i 


— 

■ Cy3 - y2) 

TAN 

■ (x2 - 

x3) 

2 


2 

L 2 J 



la 


la 


u 


U j 

r «° 


— 

■ Cy2 - yl) 

TAN 

-Cxi - 

x2) 

2 


2 

L 2 J 



L 1c 1c 

#13: Die Ankreisradien - The radii of the adj. circles 

f u fa 0 ') u f p° A u f y° 

#14: p_a := TANl I, p_b := TANl I, p_c := TAN | 

#15: [p_a := 3.864917939, p_b := 5.180911045, p_c := 11.81361004] 

"Umkreis/Circumcircle: r 

Inkreis/Inci rcle : p 

Ankreis/adj_ci rc p_a: p_a 

#16: radien Ankreis/adj_ci rc p_b: p_b 

Ankreis/adj_ci rc p_c: p_c 

p_a + p_b + p_c - p 

(p_a+p_b+p_c-p)/4: 

4 
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#17: 

ankr_a := j 

"M_ + p_a.C0S(t), 
1 

Al_ + p_a.SIN(t)j 

#18: 

ankr_b := j 

"A2_ + p_b -COS(t) , 
1 

A2_ + p_b ■ SIN(t) J 

#19: 

ankr_c := | 

"A3_ + p_c -COS(t) , 
1 

A3_ + p_c ■ SIN(t) J 

#20: 

[ankreise 

:= [ankr_a, ankr_b, 

ankr_c], adj_cnrcs := ankreise] 


Man kann die Ankreismittelpunkte auf eine andere einfache Art ermitteln: die folgenden drei Winkelsymmetra- 
len halbieren die Aufienwinkel und stehen auf die urspriinglichen Winkelsymmetralen w_a, w_fi und w_5 nor- 
mal. In deren Schnittpunkten finden sich die Mittelpunkte. 

We can find the centres of the adjacent circles in another way: the following three angle bisectors are 
perpendiculars to w_a, w_fi and w_Sand bisect the external angles. Their intersection points are the 
centres. 

#21: wn_a := [xl + |jl>(y4 - yl) , yl + pl'Gl - x4)] 

#22: wn_|3 := [x2 + |_j 2 ■ (y5 - y2) , y2 + p2'(x2 - x5)] 

#23: wn_y := [x3 + |_iB ■ CyS - y3) , y3 + p3'(x3 - x6)] 

" [wn 5,wn 15] --> A1 , [wn a,wn 15] --> A3 , [wn a,wn 5] --> A2 !" 

Mit den Koordinaten der Eckpunkte lost man die folgenden drei Gleichungen und substituiert die jeweils erhal- 
tenen Parameter. 

Using the coordinates of the vertices you have to solve the three following equations to obtain solu- 
tions for the parameters which have to be substituted into the appropriate equations. 

wn I5=wn 5 
wn a=wn 15 
wn a=wn 5 

Nun bilden wir aus alien merkwiirdigen Punkten (H, U, S, I) einen Vektor und generieren weitere ubersichtliche 
Tabellen. Im Anschluss daran berechnen wir die Gleichung der Eulerschen Geraden und zeigen, dass der 
Schwerpunkt die Strecke UH immer im Verhaltnis 1:2 teilt. 

Now we use the points H, U, S and I to create a vector and we generate more tables of results. Then 
we calculate the equation of Euler’s line and we show that the gravity centre S divides the segment 
UH in the ratio of 1:2. 



[punkte : = [h_, u_, s_, i_] , hoehen : = [h_a, h_b, h_c] ] 

[ s tree kens ymmetralen : = [ s__bc, s_ac, s_ab] , schwerlinien : = [ s_a, s_b, s_c] ] 
winkelsymmetralen : = [w a,w 15, w 5] 
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#2S: husilin := 


Hohe auf BC = a: h_a 

Hohe auf AC = b: h_b 

Hohe auf AB = c: h_c 

Str.symm von BC: s_bc 

Str.symm von AC: s_ac 

Str.symm von AB: s_ab 

Schwerl inie durch A: s_a 

Schwerlinie durch B: s_b 

Schwerlinie durch C: s_c 

Winkelsymm, von oc w_a 

Winkelsymm, von p: w_p 

Winkelsymm. von y: w_y 


#30: Die eulersche Gerade / The Euler Line 

#31: euler u_ + v- (centroid - u_) 

| u_ - centroid | 1 

#32: = 

I centroid - orthocenterl 2 


#33: Einige praktische Tabellen / Some useful tables: 


[ Dreieck [x y] 

Pkt A: [xl, yl] 

#34: angabe := , data := 

Pkt B: [x2, y2] 

L Pkt C: [x3, y3] J 

I" Winkel CAB = a: a 

Winkel ABC = p: p 

#35: winkel := , 

Winkel ACB = y: y 

L Winkelsumme: a + p + y . 

I" Seite BC = a: la 

#36: l_seiten := Seite AC = b: lb , l_side 

L Seite AB = c: 1c . 


#35 : winkel := 


Triangle [x , y] 

Point A: [xl, yl] 

Point B: [x 2 , y2] 

Point C: [x3, y3] J. 

Angle CAB = or. a 

Angle ABC = p: p 

angles := 

Angle ACB = y: y 

. Sum of angles: a + p + y 

'sideBC=a: 1 a l" 

s := side AC = b: lb 

_sideAB=c: 1 c J _ 


#37: abmessungen 


Seite BC = 

a: 

la ' 

I" side BC = a: 

la 

Seite AC = 

b: 

lb 

side AC = b: 

lb 

Seite AB = 

c: 

lc 

side AB = c: 

lc 

Urn fang : 


u 

perimeter : 

u 

Winkel oc 


a 

angle oc: 

a 




, measures := 


Winkel p: 


P 

angle $ : 

P 

Winkel y: 


¥ 

angle y: 

¥ 


Flacheninhalt f: 
Umkreisradius r: 
Inkreisradius p: 


area f: f 

radius of circ_circ r: r 

radius of incircle p: p 
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Fiir eine spezielle Angabe werden einige Ausgaben wiedergegeben. 

For one special triangle you can see some results. 

#38: [xl := -4, yl := -3, x2 := 5, y2 := 0, x3 := 0, y3 := 3] 

T Dreieck [x , y] 1 [ Triangle [x , y] 

Pkt A: [-4, -3] Point A: [-4, -3] 

#39: [angabe, data] = 

Pkt E: [5, 0] Point B: [5, 0] 

_L Pkt C: [0, 3] J L Point C: [0, 3] 

side EC = a: 5.830951894 

side AC = b: 7.211102550 

side AE = c: 9.486832980 

perimeter: 22,52888742 

angle a: 37.87498365 

#40: measures = 

angle |5: 49.39870535 

angle y: 92.72631099 

area f: 21 

radius of circ_circ r: 4,748791468 

radius of incircle p: 1,864273153 

Die Lange der Hohen - The lengths of the altitudes 

" 21.734 21-713 7.710 

#41: [lh_a, 1 h b , lh_c] = , , 

17 13 5 

#42: [ 1 h_a , lh_b, lh_c] = [7.202940575, 5.824352060, 4.427188724] 

The important points of the triangle: in German husi, in English occi: 

Important Triangle Points: [x , y] 

Orthocenter: [-0.1428571428, 3.428571428] 

CircCircCenter: [0.5714285714, -1.714285714] 

Centroid: [0.3333333333, 0] 

#43: occi = 

Incenter: [0.5651280036, 0.4868257827] 

Ad jCi rcCenter Al_: [5.464195246, 4.228712967] 

Ad jCi rcCenter A2_: [-7.324737538, 1.352913907] 

Ad jCi rcCenter A3_: [3.581128573, -12.92559551] 

The osculation points of the incircle: in German bp, in English op: 

Beriihrpunkte des Inkreises [x , y] 

an der Seite a: [1 . 5242SS627 , 2.0S5426S23] 

#44 : bp = 

an der Seite b: [-0.9860410222, 1.520938466] 

an der Seite c: [1.154662938, -1.281779020] 
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Jetzt wollen wir zeichnen - Now let us plot: 

die Seiten - the sides with [^parameters 1 

die merkwurdigen Pun He (punkte) - the important points (points) 

die Schwerlinien (schwerlinien) - the medians (medians) 

die Ankreise (ankreise) - the adjacent circles (adj_circs) with 0<pararneter<2n 

den Umkreis (umkreis) - the circumicircle (circumcircle) with 0<parameter<2jj 


#45: 

sides 

#46: 

points 

#47: 

medians 

#4S: 

adj_ci res 

#49: 

ci rc_ci rcle 



eckpunkte 

vertices 



dreieck 

triangle 



angabe 

data 



winkel 

angles 



1 seiten 

1 sides 

Seitenlangen, Lengths of sides 


abmessungen 

measures 



bp 

op 

Inkreisberuhrpunkte, osepts of incircle 


fusspunkte 

pedpoints 

Hohenfufipkte, pedalpoints of altitudes 


seiten 

sides 

Seiten, sides 

0 < p < 1 

umkreismpkt 

circcirccenter 

Umkreismittelpunkt, circumcenter 


umkreis 

circ circle 


0 < p < 2n 

hoehenschnpkt 

orthocenter 



schwerpunkt 

centroid 



inkreis 

incircle 


0 < p < 2n 

inkreismpkt 

incenter 



seitsymm 

perpbisecs 


0 < p < 1 

winkelsymm 

angbiss 


0 < p < 1 

hoehen 

altitudes 


-5 < p < 5 

schwerlinien 

medians 


0 < p < 1 

husi 

occi 

alle merkw. Punkte, all interesting points 


punkte 

points 

nur H, U, I, S, only O, C, C, I 


husilin 

occilin 


0 < p < 1 

euler 

euler 

eulersche Gerade, Euler line 


Ankreise 

adj_circs 


0 < p < 2n 

radien 

radii 
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Working sequence for the teaching unit 

The remarkable points in the triangle 

using the utility file DREIECK.MTH 


Turn on your computer and start the CAS DERIVE as usual. 

Load DREIECK.MTH as an utility file. You will not see the expressions on the 
screen but they are available for further usage. (Principle of a Black-Box!) 

Transfer Load Utility DREIECK ■ 

Read carefully the instruction how to use DREIECK.MTH (3 pages)! 

The following triangle shall be calculated and then be plotted using DERIVE : 
AABC[A(-4|-3), B(5|0), C(0|3)] 

Author [x I : = —4,yl: = —3,x2:=5,y2:=0,x3:=0,y3:=3] ■ 

Now execute the DERIVE calculations. Keep attention to the instruction (page 1 bot- 
tom and page 2 top). The expressions like angabe, abmessen, husilin have to be 

edited in the Author line without quotation marks. (That are prepared definitions 
from the Black-Box!). 

Regarding to instruction on page 2 the expressions needed for plotting the triangle 
like seiten, schwerlinien, .... ,inkreis, ankreise have to be Authored and then Simpli- 
fied or approximated. Then you can plot them. 

For the 2D-plot use an Overlay window. So the full screen is at your disposal for the 
construction. In the 2D-plot window set Options Color Auto No ■ 

Plot the sides in yellow: Options Color Plot 14 ■ 

Plot the gravity lines in light red: Options Color Plot 12 ■ 

Plot the altitudes in light green: Options Color Plot 10 ■ 

Plot the perpendicular bisectors and the circumcircle in light blue 
Plot the angle bisectors, the incircle and the adjacent circles in white! 



Important: 

Save your work as: 
MY1WORK.MTH on 
the diskette! 

Important: 

When you are ready 
with your picture then 
show it to the teacher! ! 
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An Exercise to 
Motivate Eigenvalues 

Carl Leinbach & Marvin Brubaker, USA 

In this section we will give a visual demonstration of the effect of applying a 2 by 2 matrix to a 
simple two dimensional figure. We will see this figure compressed and turned in the direction of the 
dominant eigenvector of the matrix. 

Before starting this demonstration, create a Plot window aside of the Algebra window and choose 
the Connected option for Plotting points. 

First we create the figure to which we are going to apply our matrix. Enter a 2 by n matrix that 
gives the coordinates of a house. The first row contains the x-coordinates, the second row the 
v-coordinates. 

h : = [-2 ,- 2 , - 3 , 0 , 0 , 3 , 2 , 2 , 2 ,- 2 ; - 2 , 2 , 1 , 4 , 4 , 1, 2 ,-2 ,- 2 , - 2 ] 

In order to plot the figure represented by our matrix create a vector of 2 by 2 matrices from the 2 by n 
matrix. This gives a vector of pairs of points that we can plot. NOTE a^iij is DERIVES way of de- 
noting aij. The down arrow (i) can be replaced with the phrase SUB. (Or in DERIVE 2.x versions use 
the ELEMENT(a,i , j) command). 

PMAT(b) := VECTOR([bj.lj.i ,bj.li(i+l) ;bi2ii ,bi2i(i + 1) ] ' , i , DIMENSION (b ' )-l) 

When we apply PMAT to h from above and plot the result we see: 



For the matrix, a, that we will study, generate a random matrix. DERIVE has a random number 
generator called RANDOM. RANDOM (1) generates numbers between 0 and 1. It is almost certain (but not 
guaranteed) that the following statement will generate a matrix that has a dominant eigenvalue [1] . Edit 

[RANDOM(l) , RANDOM(l) ; RANDOM(l) , RANDOM(l) ] 


(To be sure to obtain different results at any start of the file Edit and Simplify first the expression 
RANDOM(O)) . 
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#4: RANDOM® = 4054 

" RANDOM (1) RANDOM (1) " 

#5: 

. RANDOM (1) RANDOM (1) . 

Approximate this expression and then assign that result to the matrix, a. Here is one example of ap- 
proximating the above 



' 0.419901 

0.23371 " 

#6: 

. 0.665553 

0.940504 . 


Now, we assign this result to the matrix, a. We must do things in this order. Otherwise, the matrix a 
will change each time it is applied. Use the F3 key to Edit 

a := [0.419901, 0.23371; 0.665553, 0.940504] 


Let’s see what happens to the house under the action of a. 

Clear the plot window and first approximate Let’s apply a again, 

then Plot PMAT(a-h). Plotting PMAT(a A 2-h): 



Now we apply the 10-th power of a to h. Change the scale in the plot window to x = 0.02, y = 0.02 and 



It is so 'squashed' it looks like a line. If we 
apply a to h any more times we will still get 
this line. It appears to stay fixed under the ac- 
tion of a. If so, any point on the line will just 
map onto another point on the line, i.e. 

a ■ x = Sx 

for any point, x, on the line we see on the 
screen. 


Let’s find 8. Move the crosshair to the point with the x-coordinate 0.01. Our reading of this 
v-coordinate is 0.0312. 
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Now let’s see what happens to this point under 
the action of a. Author 
a . [ [0.01] , [0.0312] ] 

The result is: [0.01149076; 0.03599926], 
Dividing the x-coordinate of the above by 0.01 
we have an estimate for 8 of 
8= 1.14907. 



We can check this experimental result by finding the roots of the Characteristic Polynomial. If we 
Simplify CHARPOLY(a), we have: 


#14: CHAR POLY (a) 

2 

#15- w - 1,360404999. w + 0.2393721784 
2 

#16: N SOLVE (w - 1 . 3,60404999 -w + 0 . 2393721784 , w, Real) 

#17: w - 1.152752244 v w = 0.2076527540 


The largest root of this polynomial is w = 1.15275. 

#18: | 1.149076199 - 1.15 2752244 | = 0.003676045 


This differs from our experimental result by only 0.00368. Not bad. 
Reset the screen and define r as follows: 


#26: 



r 7T ^ 


" rr ^ 


COS 

— 

SIN 

— 



^ 6 j 


. 6 j 



( 77 

’I 

r * " 


- SIN 

COS 

— 



l 6 

J 

L 6 > 



#27: VECTOR CPMftT (r .h), k, 0, 10) 

Repeat the above exercise replacing a with r. Are the results the same? Check out the roots of the 
Characteristic Polynomial for r. 

It is very nice to use the VECTOR command to make visible the effect of ongoing multiplying the ob- 
ject matrix h by any "random" matrix a or by the "determined" matrix r. Josef 



ti] 


You can find an exciting application of dominant eigenvectors in DNL#68, page 12 
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ERF functions in DfW4 (and DERIVE 6.10) 


Another view from the 
DERIVE Window (DfW 4.0) 


Compare DfW4 and DERIVE 6 (below)! 
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Playing with DERIVE 's 
Plotting Accuracy 

Josef Bohm and Charles Reno, Euclid, USA 

When I attended the T 3 Conference in Columbus, Ohio, I met Charles Reno from Euclid, Ohio and during a 
break he challenged DERIVE to plot the implicit function 

x sin(2 y) = y cos(2x) 

Immediately we tried - in the default Accuracy 7. As we were not really satisfied with the result we increased 
the accuracy up to 9 and then, because the picture needed too much time to be built up we again decreased the 
accuracy. So we met many different patterns dependent on the plotting accuracy. It was fascinating to observe 
DERIVE ' s plotting algorithm applied on this example. 

You can follow the sequence of accuracy in the next screen shots. But one warning: if you have the intention to 
reproduce the plots then you should preserve a lot of time. 

On the next page you can find two plots made by DERIVE for WINDOWS. The plot process in DfW is very 
much faster than in the DOS version - but it seems that you don’t have any possibility to change the accuracy. 
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= File Edit Set Plot! Options Window Help 
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Solution of the problems from ACDC One 


#1: Some auxiliary functions: 

#2: [a := [0, 0], b := [1, 0], c := [1, 1], d := [0, 1]] 

#3: square(x, y, 2 , t) := [x , y, z, t, x] 


-Cx - y_ ) 

i 


#4: 

segm(x, y) := [x y] 

x- - y. 
2 

#5: 

ec line(x , y ) := y 

+ 


2 

x- - y. 
1 


rDIM(x) 

2' 

#6: 

dist(x, y) ;= J I 

(y . - x ) 


L n=l 

i i 

#7: 

1, Tsu Ch'ung Chih Approximation 

r ? l 

#8: 

*T tJ 



7 

7-x 

#9: 

ec_line(b, ej = - 



8 

8 

#10: 

Coordinates of F: 



(221 1 
#11: SOLUTIONSl (x - 1) + y = a y = ee_lme(b, e) , [x , y ] I 

' 4-413 7-^13 ' 

+ 1 

113 226 

# 12 : 

4-413 7-413 

1 

113 226 

' 1.376288347 -0 . 3292523039 ' 

#13: 

. 0.6237116526 0.3292523039 . 

#14: The first solution for x is not in [0,1]; hence we take the 2nd one: 

r 4-4i3 7-4i3 i r 4-413 r 

#15 : f := 1 , , g := 1 , 0 

L 113 226 J L 113 


#16: EXPAND(ee_lme(e, g)) 

#17: FH has the same slope 

d 

#18: — ec_line(e, g) = - ■ 


7-413-x 791- x 7 

194 776 8 


d 7-413 791 

— ec_line(e, g) = - - 

dx 194 776 


line through f with this slope: 

( 7-413 791 T 

f + ■ (x - f ) 

2 L 194 776 J 1 


7-413 791 T 7-413 


194 776 


226 8 


solve #21 = 0 tor x, giving H 


[- J , 

791 ' 

1 7-413 

h J_ J 

L L 194 

776 „ 

J 226 

h i X 

8 J 
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#24: x = 

113 


#25: h := , 0 

L H3 J 

#26: fr is approximately 3 + distance BH! 

3B5 

#27: 3 + distCb, h) = 

113 

#28 : 3 + di'st(b, h) = 3.141592920 

355 

#29 : tt 

113 

-7 


#30: 

2.667625783-10 


#31: 

1, Hobbes' Method 



2 

2 

#32: 

SOLVE (7(1 - x ) = 

1 - 7(1 - x ) , x) 


73 

73 

#33: 

X = - V X = ■ 

— 


2 

■ 73 1 ‘ 

2 

#34: 

. 22 . 



#35: Length of CE 


fd 2 "|2"| rr 

1 + — (1 - 7(1 - x )) dx = — 

ydx ) ) 6 


{ fd 2 'jZ'j 

#37: 50LVE Jl + - (1 - /(I - x )) dx = , w, Real 

- 73/2 


& 42 

#38 : w := + 

4 4 



#39: f := Lw, 1 - 40 - - w ) J 

#40: Distance between F and segment CD 

46 42 

#41: 1 - f = 

2 4 4 

#42: Coordinates of G 

g := [f , 1 - 2-(l - f )1 
#43: L 1 2 J 

f 46 43 1 76 73 3-72 5 

#44: ec_line(e, g) = - x- + - 42 - - - + + 

l 2 2 2 J 2 4 4 4 

#45: Coordinates of H 

#46: h [1, SLIBST(ec_1i' ne(e, g) , x, 1)] 

3-43 7-42 7 ' 

#47: h := 1, - Jfc + + 

4 4 4 . 






21-.J2 


D-N-L#23 


fl C D C Two 


p 45 


#48: 

tt l's app maxima 

#49: 

S-di st(h , 

c) = 

#50: 

6-cHst(h , 

c) = 

#51: 

| TT - 6-di 

st(h , 

#52: 

| tt - 6-dn 

st(h i 


9.-/3 


2-tt - 12 ■& - 9-.J3 + 21-.J2 + 9 


Another DfW Screen Dump 
(Dfw 4) 

The graphics are - like the others in this 
and earlier DSL’s - produced by students of 
Luk Van den Broeck, Belgium. 

See below the pompoen (= pumpkin) 
and the paddestoel (= mushroom) in 
DERIVE 6. 



pumpkfnCr, 5, t) := [r- SIN(s) -COS(t) , r- SIN CO ■ SINJt) , r-COS(s)- (SIN(s) + 0.5)] 
pumpki n(S) 

mush roorn(r , s, t) := [r- SIN(s) -COS(t) , r- SIN(s) ■ SIN(t) , r-COS(s) ■ (1 + COS(s)) - 5] 
mush room (8) 
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Titbits from Algebra and Number Theory 

by Johann Wiesenbauer, Vienna 

Hi folks, long time no see! Well, on second thoughts this greeting might be inappropriate in your case as 1 have 
just met quite a few of you at the second international DERIVE conference in Bonn this summer. Speaking of 
Bonn this was really a very special event! As for me it was here that I first got a lot of feedback as regards this 
column. Boy, modesty forbids me to say more than that I was really flattered! On the other hand I don’t want to 
cover up the fact that I also met people who were in the bad habit of skipping my column regularly being of the 
opinion that it is strictly for the birds. You happen to be one of them and were just about to turn over? Phew! 
Please, read on! I promise to bend over backwards to please you this time! 

Have you ever regretted the fact that DERIVE has no built-in function to solve even the simplest nonlinear sys- 
tems of equations? If you agree then the following routine might come in handy: 

SOLVE2(u,v,x,y):=ITERATE(IF(LHS(l_ SUB l)=y, APPEND(DELETE_ELEMENT(1_, 1 ),RHS( VECTOR( 

[x_,l_ SUB 1], x , LIM( SOLVE(v,x),y,RHS(l_ SUB 1 ))))), APPEND(DELETE_ELEMENT(1_,1), RHS( 

VECTOR([l_ SUB l,y_],y_,LIM(SOLVE(v,y),x,RHS(l_ SUB 1))))), 1_), l_,ITERATE(SOLVE(g_,IF( 
PRODUCT (v_-x,v_, VARlABLES(g_))=0,x,y,y)), g_, ITERATE(RHS(v)*QUOTIENT(u_,LHS(v))+ 
REMAINDER(u_, LHS(v)),u_,LHS(u)-RHS(u)),l)) 

Here u and v are two polynomial equations in x and y where v should be the “simpler” of the two equations. It is 
supposed to work whenever 

RED(u,v):=ITERATE(RHS(v)*QUOTIENT(u_,LHS(v)) + REMAINDERS, LHS(v)),u_,LHS(u)-RHS(u)) 

(cf. DNL #20) yields a solvable expression in one variable only. 

Using these functions example 1 from DNL #21, p.39, would now look like this: 

RED (x A 2-x*y-y=113, y A 2-x*y+x=-41) =126* (y+13) / (y-1) A 2-70 
SOLVE2 (x A 2-x*y-y=113, y A 2-x*y+x=-41) = [ [15,7] , [-61/5,-16/5] ] 

I used what I call the iterate(...,l)-trick to “combine RED and solve”, as Josef put it in the article quoted above, 
and it might be of use in other cases as well: 

SOLVERED(u,v):=ITERATE(SOLVE(u_,VARlABLES(u_)),u_,RED(u,v),l) 

Actually, it was this preliminary version of SOLVE2 which I distributed to some people at the conference in 
Bonn. It is still useful in some cases like in my example 2 in the DNL #21, p.39: 

P(v) :=(a*(b-v)+r*t*v A 2)/ (v A 2* (v-b) ) 

SOLVERED (DIF (P (v) , v) , DIF ( P (v) , v, 2 ) ) = [v=3*b, v=inf , v=-inf ] 

Here SOLVE2 won’t work because of a very nasty bug in the current version 3.14e of DERIVE. Which one? 
Well, DERIVE commutes sometimes the operations DIF and LIM where this isn’t allowed, e.g. in 

LIM(DIF(xy,x),y,x)=2x 

There are other cases where SOLVE2 should work in theory but doesn’t in practice because of a memory over- 
flow e.g. 

SOLVE2(x A 2-2*x*y-3*y=l,y A 2-3*x*y+2*x=-4) 

In these cases the modified routine 

NSOLVE2(u,v,x,y):=ITERATE(IF(LHS(l_ SUB l)=y,APPEND( DELETE ELEMENT (1_, 1 ),RHS( VECTOR( 
APPROX([x_,l_ SUB 1]), x_, LIM(SOLVE(v,x),y,APPROX(RHS(l_ SUB 1)))))),APPEND( 
DELETE_ELEMENT(1_, 1 ),RHS( VECTORf APPROXf [1_ SUB l,yj), y_,LIM(SOLVE(v,y),x,APPROX(RHS( 
1_ SUB 1)))))),1_),1_, ITERATE(SOLVE(g_,IF(PRODUCT(v_-x,v_,VARIABLES(g_))=0,x,y,y)), g_, 
ITERATE(RHS(v)*QUOTIENT(u_,LHS(v))+REMAINDER(u_,LHS(v)),u_,LHS(u)-RHS(u)),l)) 

should give you at least approximate solutions (n stands for numeric!). Unlike the usual numeric procedures you 
will get all the solutions (even complex ones if any!): 

NSOLVE2 (x A 2-2*x*y-3*y=l , y A 2-3*x*y+2*x=-4 ) = 

[ [3.58397, 1 . 1649] , [-1 . 68707, -4 . 93439] , [-1 . 04 845-0 . 353058*#i, -0 . 415254 + 

0 . 495038*#i] , [-1 . 04845+0 . 353058*#i, -0 . 415254-0 . 4 95038*#i] ] 

(440 sec. on a 486DX2/66; 575 sec. on a 486/33; 110 sec. on a Pentium) 
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You might have noticed that I have always omitted the third and fourth parameter so far. In fact, this is allowed 
if your variables are x and y which is recommended. If you use other variables than x and y you should make 
sure that they come before any constants in the order of variables (if necessary use the command “Manage Or- 
der”) e.g. 

SOLVE2 (a A 2+b A 2=l,b-r/s*a=-l,a,b) = 

[ [0,-1] , [2*r*s/ (r A 2+s A 2) , (r A 2-s A 2) / (r A 2+s A 2) ] ] 

It is not easy to give general rules when SOLVE2(u,v,x,y) is applicable and therefore it is up to 
you to experiment a little. Even so, the following two cases deserve to be mentioned: 

• SOLVE2(p(x,y),ax+by=c), where p(x,y) is a polynomial of degree at most 4 and a,b,c are real numbers. (In 
the case a=0 reverse the order of x,y by means of “Manage Order”!) In particular you could use 
SOLVE2(p(x,y),x=c) or SOLVE2(p(x,y),y=c,y,x), respectively, to calculate points on the curve which is im- 
plicitly given by p(x,y)=0. 

• When dealing with two conic sections try to arrive at SOLVE2(ax A 2+bxy+cy A 2+dx+ey+f,gy A 2+hx+iy+j) 
with h^O or alternatively at SOLVE2(ax A 2+bxy+cy A 2+dx+ey+f,gx A 2+hx+iy+j,y,x) with i^O (again reverse 
the order of x,y in the latter case!) In “school examples” this special form of v can be virtually always 
achieved by forming suitable linear expressions of the original equations u and v. Since SOLVE2 might take 
quite long in these cases (or even lead to a memory overflow) you should always try SOLVERED or 
NSOLVE2 first! 

Let me use the remaining space to deal with another “most missing” feature of DERIVE which is sorting a given 
list of real numbers. Those lists occur in mathematics on many occasions. One example that immediately comes 
to my mind are the Farey numbers of order n for a given natural number n, i.e. the positive fractions not greater 
than 1 with a denominator at most n ordered by their size. Assuming that a routine SORT(v) is already available 
the definition of FAREY(n) could look like this: 

FAREY(n)-SORT((ITERATE(IF(a_=b_ AND b_=n,[a_,b_,c_], [IF(a_=b_,l,a_+l), 

IF (a_=b_,b_+ 1 ,b_),IF(GCD (a_,b_)= 1 , APPEND(c_,[a_/b_]),c_)]),[a_,b_,cJ,[l,l,[]])) SUB 3) 

As for SORT we could take in this case 

SORT(v):=REVERSE_VECTOR(LIM(TERMS(v*VECTOR(x A k_,k_,v)),x,l)) 

Actually this is nothing else but the sort-routine from my “Titbits(7)” which I believed only sorted different 
nonnegative integers but which can be taken for different rational numbers as well! Nevertheless I will keep my 
promise and give you the definition of a sort-routine that sorts arbitrary real numbers: 

SORT(v):=REVERSE_VECTOR(LIM(TERMS(EXPAND(VECTOR(k_*x_ A APPROX(k_),k_,v)* VECTOR) 
y_ A k_, k_,D IMEN SION(v)))), [x_,y_],[l,l])) 

(It goes without saying that in both cases the corresponding definition for DSORT(v), which sorts in descending 
order, is the same but without applying REVERSEVECTOR.) It works as the following example shows: 

SORT ( [pi, SQRT (3) , SQRT (2) , #e, SQRT (5) , SQRT (2) , #e] ) = 

[ SQRT (2) , SQRT (2) , SQRT (3) , SQRT (5) , #e, #e,pi] 

I hear you saying that there has already been published such a function in the DNL #22,p.40, by Sergey 
Biryukov. Ah, I would have liked to dodge that issue because it forces me to admit that I felt totally at a loss 
when I studied Sergey’s solution. But remembering a story dealing with Alexander the Great when he was in a 
quite similar situation I cut Sergey’s implementation into two halves throwing away one of them and making 
some changes to the other. Here is the outcome of my doing: 

SSORT_AUX(v,n,p):=VECTOR(IF(MOD(i+p,2)=0,IF(i=l OR v SUB (i-l)<= 
v SUB i,v SUBi,v SUB (i-l)),IF(i=n OR v SUB i<=v SUB (i+1), v SUB i,v SUB (i+l))),i,n) 

SSORT(v):=ITERATE(SSORT_AUX(SSORT_AUX(w,DIMENSION(v),0),DIMENSION(v),l),w,v) 

Just one small example to give you an impression of the calculation times: For n=50 there are 774 Farey frac- 
tions to sort. It took 7 sec to calculate them and 1.1 to sort them with the first sort-routine, 1.5 sec with the sec- 
ond one and 901.4 sec by means of SSORT(v). 
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I should also note that I am currently using a Pentium PC with 166 MHz and 32 MB Ram. It’s really a powerful 
machine! 1 used it to calculate the number 

parts(1840926)=245...373 (1505 digits!) 

in 12 h 43 min which has been proven to be prime by Morain quite recently . It took DERIVE another 394.3 sec 
to confirm this. See you next time!(?) 

On the diskette Johann gave the solution of the example from DNL#21 using his SOLVE2 routine: 

2 2 2 

#19: SOLVE 2 (a +b +c =9, a + b + c = 3, a, b) = 

V3 • V (c - 3 ) • V ( -c - 1) + c - 3 V(c - 3) • (V3 • V (-c - 1) - V(c - 3)) -| 

2 2 
V3 • V ( c - 3) • V (-c - 1 ) - c + 3 V3 • V ( c - 3) • V (-c - 1 ) + c - 3 

2 2 
V3 • V ( c - 3) - V (-c - 1 ) + c - 3 

#20: a := - User 

2 

V (c - 3 ) • ( V3 - V ( -c - 1) - V (c - 3)) 

#21: b := User 

2 


3 3 3 3 2 

#22: a + b + c - 24 = 3-(c - 3-c + 1) 

3 3 3 

#23: SOLVE (a + b + c = 24, c) = [c = 1 


User=Simp (User) 


1 


2 - COS 



#24: "The following equation together with a A 3+2b A 3=24 shows that" User 

#25: "no two variables, like e.g. b and c, may have the same value:" User 

2 2 3 0 

#26: SOLVE 2 (a + 2-b = 9, a + 2 -b = 3, a, b) = [] ] User=Simp (User ) 


#27: "Therefore apart from permutations this is the unique solution:" User 


#28: 

#29: 

#30: 

#31: 

#32: 

#33: 

#34: 


|"a := 1 - 2 • SIN f — — ] , b : 

= 2 

•COS 

r— 1 

\ — 1 
+ 

c := 1 - 2 • COS [ 

- 2 • n - 1-1 

User 

L L !8 J 

"Let's check all the terms 

in 

L 9 J 

question : 

II 

l 

- 9 4 4 

User 

Trigonometry := Collect 

r 2 2 2 

L a + b + c, a + b + c , 

3 

3 

3 

4 

4 4-i 

+ b + c J 

User 

a 

+ b 

+ C , 

a 

User 

[3, 4 • COS f— — — 1 - 2 • V3 • COS f 

n 

-1 - 

2 -SIN 

[- n 

-1 + 9, 24, 69] 

Simp (#31 ) 

L L 9 J L 

"DERIVE can't cope with a A 

■ 18 
2 +b 

A 2 +c x 

'2 but 

L 18 
by 

J J 

approximating we 

> arrive at" 

[3, 9, 24, 69] 






Approx (#31) 
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DIRA - A TI-92 Program for Discussing Rational Functions 

by W. Propper, Nurnberg, Germany 

When using the TI-92 more often to do the necessary computations for discussing functions one soon 
sees what a boring activity this is. The same operations have to be done again and again: 

Enter the term of the function, store it as f(x), let the derivatives be computed and store them with 
other names, compute zeros of the function and its derivatives and so on. At the end you can get a nice 
plot. 

What was more obvious than using the programmability of the TI-92 to pass all the routine works to 
the calculator. I was inspired to do so by a little contribution by J. Bohm in DNL #15 from Sept. 1994. 
There a Derive-utility was introduced, in which only the function’s term had to be entered and Derive 
did the whole discussion. 

This contribution shows in a comprised form the possibilities of DIRA when discussing rational func- 
tions. 


After inputting diraQ ( ENTER) on the home 
screen of the TI-92 you see this opening screen: 



Discussion of Rational Functions 

[FI 3 tils Help and instructions 
[ F 2 3 5 Ente r f u n c t i o n 

< o t h e r r-i e n us inaccessible) 

fig ■ i 

Looking into the menus an expert foresees 
which mathematical activities take place in each 
item. 

Algebraic investigations with ( F3) : 



fig ■ 4 


Tools are found in (FI): 


' fit 

Tools 


2: fll INew 

3: Fluto/Flpprox Mode Tonal Functions 
4: End 

5= About DIRfl 

TFTTTTT! Help and instructions 

fig ■ 5 


To be able working with DIRA a function’s 
term must be entered first via ( F2) : 


fig ■ 2 

Investigations by means of calculus with ( F4) : 


Enter numerator and denom. as poly- 
nomials of n ot higher than 5th d egree. 

Numerator: I_T2^3+Bx22jj2x I 

Denominator: 19 ( x- 1 !> ^21 1 

(The input lines can be edited. ) 


ESC=CHNCEL 




fig- 7 


Graphs and tables with ( F5) : 
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By ( F3) ( 6) all algebra-items are run through: Extrema and inflection points can be found in 



fig- 9 

Finally the graphs of the function and its first derivative (dotted): ( F5) ( 2) and ( 3) : 



The program with a more detailed documentation will be found on the annual DNL-diskette. By 
means of the TI-Graph Link it can be transferred from the PC to the TI-92. From autumn 1996 you 
can find it in its most recent version on the home page of my school with the web address 
http://www.nuernberg.de/ver/schg/nk/ . 

Many thanks Wolfgang for your wonderful "Christmas gift". Maybe I’ll find a way to put your files in a 
read- and printable form on the diskette, too. Then our friends who have not GRAPH LINK available 
can also benefit from your program. See the first lines of DIRA 's HELP subroutine: 

0 

P r g m 
Cl rl 0 

Di sp " Hel pscreen ( 1) " 

Di sp " " 

Disp "DIRA is a program which does the routine-" 

Disp "work when Discussing RAtional functions." 

Disp " " 

Disp "In menu [F3] algebraical methods like" 

Disp "computing domains, zeros, poles and" 

Disp "asymptotes are applied." 
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Determining Equations of Conic Sections 

by M. Brubaker & C. Leinbach, USA 


Introduction: 

Given the following figure, find the equation of the circle. Where is its center located? 
What is its radius? 


Tl^wr - f:= fh 

t j— Zoom Tr ace ReGraph 

' FET Y FfiT Yff a Y 1 

Math Draw t y 

( 

) 



1 MAIM RAD HUTD FDL I 


The strategy is obvious. Use the cursor to digitize some points and use these points to de- 
termine the equation of the circle. How many points do we need to digitize? The answer be- 
comes clear if we look at the standard equation for a circle. 

(. x-h ) 2 + (. x-k ) 2 = r 2 

There are three unknowns, the x-coordinate of the center, the y-coordinate of the center, and 
the radius. We can find the equation by obtaining three points that lie on the circle, substitut- 
ing, and solving the resulting three equations for the unknowns, h, k, and r. This process is 
straight forward, but, because of the fact that the unknowns appear to the first and second 
power in the resulting equations, tedious. We propose a different strategy that can be easily 
extended to find the equations of other conics in standard form. 

Let’s suppose that (x ,y) is a point on the circle, then 

{x-h) 2 + {x-k) 2 = r 2 

which we can expand to 

x 2 + v 2 - 2hx-2ky + {k 2 + h 2 -r 2 ) = 0. 

This has the general form of the homogeneous equation 

A{x 2 + v 2 ) + Bx + Cy + D = 0. (1) 

Because we are dealing with the equation of a circle, x 2 and y 2 terms have the same coeffi- 
cients. Thus, when we substitute for x and y we have a linear system with unknowns A, B, C 
and D. Have we increased the number of points needed to find the equation? The answer is 
no. The symbolic capabilities of our CAS will show that this the case. 

Let’s get specific. The figure depicted at the beginning of our discussion was chosen 
so that we could easily identify that the points (-1,1), (6,8) and (7,7) lie on the circle. Substi- 
tuting these points into equation (1) above we have: 







(2) 
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2A - B + C + D = 0 
100 A + 6B + 8C + D = 0 
49A + IB + 1C + D = 0 


This still leaves us short one equation, or does it? We still have equation (1) and our CAS. If 
we think of (x,y) as a general point on the circle, we have a system of four homogeneous 
equations in four unknowns. From the theory of equations we know that this system will have 
a nontrivial solution if and only if the determinant of the coefficient matrix is zero. We enter 
this matrix using the Matrix editor of the TI-92. 



F2 

Plot Setup 

' F3 ' 

Cell 


■ f 

' FSt ' 

Util 

' v? ' 
Stat 


MHT 

HxH 






cl 

c2 

c3 

c4 

c5 

l 

m am 

X 

y 

1 


2 

2 

-1 

l 

1 


3 

9S 

7 

7 

1 


4 

100 

6 

S 

1 


5 






6 
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rlcl=x A 2+v A 2 

MAIN DEG flUTD POL 


In the following figure we take the determinant of this matrix (which we called conic). Setting 
this result equal to 0 yields the equation of the circle. Completing the square will reveal cen- 
ter and radius. 


To graph the circle, it is best to use polar equations. Once again the CAS makes this 
easy. We use local substitution with x = r.cos(0) and y = r.sin(O) and solve the resulting ex- 
pression for r (*0). This generates the graph at the start of this article. 


rFiT’mi 

f F2T 

Algebra 

Calc 

Other 

f rs 1 
PrgrolO 

Ffi 

Clear a-z... 


■ deU conic) 14 [x 2 - 6 x + y - (y - S)j 

■ 14 (x 2 - 6 ■ x + y (y - S)) | x = r ■ cos(&'} and ► 


14 ■ r (r - 2 (3 cost! + 4-sin(0))) 

■ solve( 14 ■ r (r - 2 (3 ■ cos(O) + 4 ■ sin(O))) = ► 

r = 2-(3-cos(0) + 4 sin(0)) or r = Q 

■ 2-(3-cos<e) + 4 sin(e)) * rlO) Done 

Graph 2*<3*cosO>+4*sinO>> 

MAIN RAD AUTD PPL lfl/30 

This same idea can be extended to the other conics. The general form for a parabola is 

Ax + B x + C y + D — 0 or Ay -\-By-\-Cx-\-D = 0 

and for ellipse or hyperbola (depending on the location of the points) is 

Ax + B y + C x + D y + E = 0 

Locating three points on the graph of a parabola or four on the graph of an ellipse or hyper- 
bola and entering the appropriate matrix will quickly generate the equation of the conic in 
question. 

Locating three points on the graph of a parabola or four on the graph of an ellipse or hyper- 
bola and entering the appropriate matrix will quickly generate the equation of the conic in 
question. 







